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I. INTRODUCTION 


Electron probe micro analysis provides a rapid 
method for the chemical analysis of volumes of the 
order of 10“Hcm3 (lOum^) ; however, the process of 
converting the data, which come from the probe in 
the form of measured x-ray intensities, into ab- 
solute concentrations is long and cumbersome when 
carried out by conventional methods. Since the 
task is obviously one for a computer, /Fortran-63 
programs have been written that not only perform 
the necessary absorption and fluorescence correc- 
tions but also carry out almost the entire data 
handling operation. With these programs, the 
amount of data that would come from a normal eight 
hours of microprobe operation can be processed in 
a little more than one minute of CDC 3600 computer 
time, at a fraction of the cost that would be re- 
quired for an individual to do the work using a 
desk calculator. 


In writing these programs we have tried to 
satisfy three basic aims: A) versatility, B) use 
of the most appropriate correction procedures, 
and C) elimination of all hand calculating and 
data sorting. 


A. Versatility 

The programs are suitable for any type of prob- 
lem from the analysis of a binary alloy to that of 
a complex mineral system involving as many as 
twenty elements. They accept any reasonable amount 
of data in any sequence in which the probe operator 
wishes to gather them and furnish the results in 
graphical as well as printed form. 

Part of the flexibility of this method comes 
from the fact that we have written not one but two 
programs, which may be run together or separately. 
The first program, EMX, takes raw data from the 
probe, corrects for deadtime, drift, and back- 
ground, and calculates the ratio of observed in- 
tensities (sample over standard) times the concen- 
tration of the element in its standard; the results 
from EMX are referred to herein as the first ap- 
proximation? The output from Program EMX can also 
be written on magnetic tape or punched on cards for 
use as input to Program EMX2, 

The second program, EMX2, makes corrections 
for absorption and fluorescence by characteristic 
radiation. Additional corrections, such as that 
for atomic number effect, can be added later when 
a satisfactory correction procedure is determined. 


Note that our first approximation is not the 
"premiere approximation" of Castaing (2) but re- 
fers only to the output of the first program (EMX). 
The first approximation is thus the concentration 
of the particular element in its standard times 
the ratio (sample to standard) of observed inten- 
sities corrected for deadtime, drift and back- 
ground but not for absorption or fluorescence. 


As many as 800 probe data points (each point 
being *made up of from one to 100 separate read- 
ings) may be included in one problem, and any 
number of problems may be grouped in one computer 
run. Within each problem, the probe operator may 
analyze up to nine different elements in the 
first program. For the second program, output 
from two or more problems may be grouped together; 
this allows the analysis of samples containing as 
many as 20 elements. 

Some savings in time and computer costs are 
made if the two programs are run together, but 
often this is not desirable. The analyst may 
wish to look at the first approximations before 
deciding if the results are worth correcting 
further. Or he may have analyzed only some of 
the elements present In the sample and wish to 
measure or estimate the concentrations of others 
before correcting for absorption and fluorescence. 
Another practical advantage of the two-program 
method is that it effectively increases the 
storage capacity of the computer and enables one 
to handle a large amount of data. 

Program EMX2 can be modified slightly to run 
"in reverse." If two cards in the main program 
are changed** true concentrations can be entered 
as input data; the output then shows what the 
observed intensity ratios will be if the correc- 
tion procedures used are both valid and complete. 

B, Appropriate Correction Procedures 

When one is writing a computer program and 
changing from manual to automatic data processing, 
it is an opportune time to re-examine all the tech- 
niques and formulas used in the correction of 
microprobe data for quantitative analysis. We 
would be failing to make use of all the advantages 
offered by electronic computers if we continued to 
approach problems from the point of view of a man 
with a slide rule. Several approximations that 
were previously tolerated because more rigorous 
calculations would have been too time consuming 
are no longer justified. 


For purposes of this discussion, a "reading" is 
defined as one line of data as it comes from the 
probe and includes the counting time, sample cur- 
rent, counts on three spectrometers, a six-digit 
code, and, if desired, an identifying label 
(NAME) . A "data point" is made up of one or more 
readings of the same sample. All the readings in 
a data point will be added together to find the 
average counts per second for that data point. A 
"data set" includes consecutive probe data points 
taken without a break. A background or standard 
peak data point may pertain to all the sample 
data points within the set, but it will not be 
used to correct data in another set. Therefore, 
each data set must contain the necessary standard 
and background readings. 

*Change Statement 11 to OCON(I) = 0C0N(I)/(SCF(I) * 
SCON(I)). Change fifth statement after Statement 
400 to C0N2 = OCON(I) *(FCHI*(1.0 + FL)). Run with 
IMAX = 1. 
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Consider, for example, the "simple" formula 
for fluorescence corrections derived by Wittry 

CD: 
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This formula was ingeniously arranged to take the 
form of three separate factors, each dependent 
upon only one variable. This is a considerable 
advantage to a person making corrections by hand. 

A computer cannot readily read values from a 
graph, but it can evaluate complicated functions 
extremely quickly. We have worked out a formula 
for fluorescence corrections which combines 
Wittry's voltage dependence factor with the 
formula of Castaing (2) , which takes into account 
the effects of depth of penetration in the sample. 
Theoretically the combination formula appears more 
satisfactory than the others, but whether it is a 
significant improvement remains to be seen. 

For the absorption correction the computer 
program now uses Philibert's (3) formula with 
Duncumb and Shields' (4) effective sigma, but it 
is hoped that a correction method based on sound- 
er theoretical considerations can be developed. 

One of the major advantages of computer handling 
of microprobe data is that it removes much of the 
drudgery from the process of accumulating large 
quantities of analytical data on known samples 
with a view to evaluating correction methods. The 
effects of making changes in the correction formu- 
lae can readily be seen by simple recycling of 
the data through the computer. Evaluations of 
the modified fluorescence correction used in EMX2 
and of the absorption correction are being attempt- 
ed by the analysis of simple known phases, particu- 
larly binary alloys. The ability to run the pro- 
gram "in reverse" and generate observed intensities 
from known compositions proves very valuable in this 
project. 


C. Elimination of Manual Data Processing 

Our final aim, which was actually the most dif- 
ficult to carry out, was to eliminate all hand cal- 
culation as well as the need to transfer data to 
some special form. On the other hand, we did not 
wish to restrict the probe analyst to following 
any certain sequence of collecting data. For 
Program EMX, output from the probe can be on punch- 
ed paper tape, punched cards, magnetic tape, or 
typewritten pages which can go directly to the key- 
punch operator. The data do not need to be re- 
arranged in any way — assuming, of course, that the 
probe operator did not forget to include the nec- 
essary standard and background readings or make 
some other error. The program user need prepare 
only one data card for each element analyzed in 
the problem (the Element Card, see Appendix 2) 
plus one card giving general information about 
the entire problem (the Identification Card). 

The technique by which this has been accom- 
plished is really the key to the entire EMX pro- 
gram. It consists in the assignment of a six- 


digit code number to each line-of probe output at 
the time it comes from the probe. Two digits 
of the code refer to each of the three spectrom- 
eters, one indicating which element of the sample 
is being analyzed and the other specifying the 
type of reading — standard peak, sample peak, 
standard background, sample background, or dead- 
time data (see Sec. II-A) . These code numbers 
could be written by hand on the data sheet that 
comes from the probe, but we have actually modi- 
fied our microprobe by adding a "green box", 
which is located conveniently for the probe 
operator and connected to the printout device. 

By setting the six switches on this box, the 
operator presets the numbers in the code column 
before the data is printed or punched. 

Because each line of data from the probe is 
identified by this code number. Program EMX can 
accept data in any order in which it is run, with 
no need for it to be arranged into a particular 
sequence. The operator can take readings on stan- 
dard peaks first, last, and/or anywhere in the 
middle of the data. Background data points can 
appear before and/or after each peak data point, 
or one background data point (or pair of points) 
can be subtracted from all peaks for the same 
element. The user may take one or more readings 
on each point; in the latter case the computer 
averages the readings for each data point and 
computes the standard deviations. It is unneces- 
sary to take the same number of readings on each 
type of data; one could, for example, take two 
readings on each sample peak, a large number of 
readings on the first standard peak ( to check 
probe stability) , and some other number of read- 
ings on later standard peaks. Data points on 
samples in any one problem, however, must contain 
the same number of readings. 

Further routine paper work has been eliminated 
by permanent storage within Program EMX2 of all 
the various constants necessary for carrying out 
corrections involving any elements. Of course, 
this data is not really permanent; it can be 
changed easily when better values become known. 

The program does not store all the values of u/p, 
o, K- fluorescence yield, critical excitation 
potential, and absorption jump ratio, but rather 
the constants necessary for calculating them. 
HeinrichfeO) method of calculating mass absorption 
coefficients and absorption jump ratios is used, 
so his values for absorption edges plus the nec- 
essary coefficients and exponents are stored. Also 
included are the atomic weights and characteristic 
wavelengths of all the elements (see Sec. Ill- A) . 

We chose this method because our computer does 
not have enough quick access storage space to hold 
vast numerical tables, and it can re-evaluate 
functions more economically than it could look up 
the answers on a tape. For a typical problem in- 
volving six elements plus oxygen, the CDC-3600 
computer took 10 milliseconds to decide which ele- 
ments will fluoresce which other elements. Calcu- 
lation of all the required physical constants, 62 
separate variables in all, took 38 milliseconds. 
This procedure is obviously faster and more eco- 
nomical than looking up numbers in tables and re- 
moves the possibility of human error in so doing. 


* * * 
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The remainder of this report is devoted to a 
detailed description of the two computer programs 
and how to use them. Neither the history of 
microprobes nor the theory behind methods of cor- 
rection for absorption and fluorescence effects 
will be discussed. For a recent, comprehensive 
treatment of the latter topic, the reader is re- 
ferred to Philibert (3). An English translation 
of this paper is presently being prepared by our 
group and will be available soon. 

Appendix 1 contains Fortran listings of the 
two programs, and data input specifications are 
in Appendix 2. Appendix 3 is made up of flow 
diagrams of most of the routines in the programs 
for those who wish to understand the logic of the 
programs and/or modify them. 

Programs EMX and EMX2 are written in CDC 
Fortran-63 and should work in their present form 
on most CDC 3600 or 1604 computers. After con- 
version to IBM Fortran they would be suitable for 
use with other large computers such as IBM 7094.* 

*Program EMX requires about 23,000 storage loca- 
tions, of which 11,000 are for the program itself 
and 12,000 for the data required for 800 data 
points. Program EMX2 requires 12,650 locations— 
2,000 for the program, 5,000 for physical constant 
data, and about 12,000 for data involved in a prob- 
lem of up to 500 samples with up to 20 elements 
each. The programs were written without any great 
effort to save space, and storage requirements 
could be greatly reduced by modified programming 
and reduction of the amount of data handled at 
once. However, modification of the programs for 
a really small computer would be extremely dif- 
ficult. 


Such conversion is relatively simple and can be 
done by any experienced programmer.** 

The programs were designed for use with an 
Applied Research Laboratories electron probe 
microanalyser. Because the data input routines 
handle data only in the format in which it comes 
from our microprobe, further modifications will 
be required by some users. The graph plotting 
routines, which were developed by California 
Computer Products, Inc., for use with their 
Calcomp 165 incremental plotter, will also have 
to be replaced with different routines if the 
program is to be used with another type of equip- 
ment. However, the basic method of data proces- 
sing embodied in programs EMX and EMX2 should be 
applicable to a wide variety of situations, and 
it is hoped that they will be useful to most 
microprobe analysts who wish to use computers 
for data reduction. 


**An IBM Fortran-IV version of these programs is 
presently being developed by the General Atomic 
Division of General Dynamics, P.0. Box 608, San 
Diego. For a listing, contact Tom West. 


II. PROGRAM EMX 


A. Data Coding and Input 

Data from the electron mi crop robe is read into 
the computer from either magnetic tape or punched 
cards. Each card (or tape record) contains one 
line of output from the probe and consists of 36 
digits plus up to eight alphanumeric characters in 
the format described in Appendix 2. Subroutine 
EMXINPUT must be modified if data is to be in a 
different format. 

The six- digit K0DE* which is part of each line 
of probe output is the key to this method of han- 
dling data in almost any sequence in which they 
happen to be run. The first, third, and fifth 
digits of the KODE specify the type of reading on 
the first, second, and third spectrometers, respec- 
tively, as follows: 

0 = sample peak 

1 - sample background 

2 = standard peak 

3 = standard background 

4 = deadtime data 


*In this discussion, names of variables given in 
capital letters are the actual names used in the 

Fortran program. 


The second, fourth, and sixth digits are the num- 
bers assigned to the elements being considered in 
that reading. An element number of zero would 
indicate that the reading for that spectrometer 
should be ignored. 

A special Break Card, with all zeros in the 
KODE column, signifies the end of a data set. 

This card is inserted in the data when there is 
a time lapse between groups of readings in the 
same problem. Another special card, with all 
nines in the KODE column, signifies the end of 
all data for a problem. 

Data can come from the probe directly on punch- 
ed cards, magnetic tape, punched paper tape, or 
typewritten data sheets. The best method for any 
particular laboratory would depend upon the type 
of operation and what facilities are available. We 
find that m most cases it is preferable to have 
punched cards (which are punched from our type- 
written data) so that one can easily eliminate 
faulty data and correct operator errors. For 
occasional long routine projects, though, we use 
punched paper tape which is converted directly to 
magnetic tape. 
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Figs. 1 and 2 are examples of typewritten data 
sheets as they come from our microprobe. They il- 
lustrate two different types of problems: 1) anal- 
ysis of six elements of a mineral and 2) analysis 
for a trace element preceded by a deadtime deter- 
mination. The comments explain the use of the 
KODE and possible methods of grouping data. 

For oxide and oxy anion systems the amount of 
oxygen associated with the various elements is cal- 
culated by the computer, assuming a known specific 
stoichiometric relationship and in Program EMX2 
these are modified as the element concentrations 
are modified by the correction procedure. The 
assignment of oxygen to a particular element is 
a function of its valence state and this factor 
(see Appendix 2) is decided by the operator and 
fed into the program. The program can also handle 
a second element (for example silicon in silicates) 
in like manner to oxygen.* In pyroxene analyses, 
for example. Mg, Fe and Ca can be analysed simul- 
taneously and the corresponding silicon and oxy- 
gen contents for each data point can be calculated 
by the computer. These components which are not 
experimentally measured, but determined from 
assumed stoichiometric relationships, are in the 
following referred to as "stoichiometrically deter- 
mined components" , abbreviated as s.d. components. 


The TIME value output of our microprobe is the 
counting duration of each reading, but for the 
drift correction one really needs to know the time 
that has elapsed between the beginning of the data 
set and each data point. (It must be remembered 
that the time scale cannot continue over more than 
one data set because the Break Card between sets 
may indicate a time lapse of hours.) Therefore 
each data set is considered to begin at TIME = 0, 
and the TIME of each data point is changed to the 
number of seconds elapsed between the beginning 
of the data set and the median time of that point. 

1. Deadtime Correction' 


The equation used for the dead time correction 
is 


COUNTS = 


PC 

1.0 - OC-TAU 


where 0C = observed counts and TAU = dead time in 
seconds. 

A value of TAU may apply only to one of the 
channels and may vary with time. Thus TAU values 
are not built into the program but must be read 
into the computer with data for each problem. 


B. Data Processing Procedures 

Data input is handled by Subroutine EMXINPUT. 
First a single line of data is read, and the com- 
puter uses the KODE to determine the type of data 
point that the reading is part of. Then any re- 
maining readings in that data point are read. At 
this stage the computer checks to see that data 
are grouped according to the specifications on the 
Identification Card (see Appendix 2) for the prob- 
lem and that dimension limits are not exceeded. If 
errors have occurred, the data cannot be processed 
correctly; so the computer skips over the data for 
that problem and goes on to the next problem, if 
any. 

If the data point was made up of more than one 
reading. Subroutine EMXSUM is called to average 
the readings and compute the standard deviation by 
two methods: the root mean square (r.m.s.) and the 
square root of the number of counts (Vn) . If the 
data point is a standard peak, it is important that 
the readings do not scatter over too wide a range 
because a large error in measured intensity of the 
standard would make all the analyses in the prob- 
lem worthless. But when automatic data processing 
is used, the data may go from probe to computer 
without inspection. Thus the computer checks stan- 
dard deviations on all standard peaks to see that 
the r.m.s. value does not exceed the /N value by 
more than a factor of SLIMIT, The SLIMIT factor 
is input with the data for each problem, and we 
suggest a value of 2.0 to 3.0. Setting the limit 
too low will cause the rejection of some good data, 
and setting it too high will have the effect of 
removing this check on data scatter. 

For data points other than standard peak, the 
larger of the standard deviation values calculated 
is retained, and COUNTS are converted to counts 
per second. When all data for one problem have 
been read, control returns to the main program. 

— — 

If s.d. components are elements other than oxygen 
and silicon, subroutine OXYGEN must be modified. 


If TAU values are not known, they may be com- 
puted from experimental data by Program EMX, 
either as a separate problem or as part of any 
other problem. Dead time data, identified by the 
numeral "4" in the KODE column, consists of a set 
of readings on one particular substance over a 
wide range of sample currents. For each data 
point, counts per second are divided by sample 
current, and the resulting values are fit linear- 
ly by a least-squares method as a function of 
sample current, i s . This is the equivalent of 
fitting the counts-per-second values with a para- 
bola which goes through the origin, described by 
the formula Ai s + Bi'| = 0. Then TAU = -B/A^ sec. 
(Figs. 3 and 4 illustrate output from a deadtime 
determination) . 

Since dead time is a function of photon energy 
and consequently of the elements analyzed, our 
first through t in writing Program EMX was to pro- 
vide a simple means of recalculating TAU at fre- 
quent intervals. We soon realized, however, that 
to obtain a reliable value for TAU one must take 
many readings at many different sample currents 
in order to obtain adequate statistics. The cur- 
rent must be measured with an extremely sensitive 
and properly calibrated meter. Thus we have found 
that an occasional but thorough determination of 
dead time is preferable to frequent quick determina- 
tions. 

Each Element Card (see Appendix 2) of the in- 
put data for Program EMX must contain information 
about the dead time correction to be made. TAU, 
if known, is entered in E format; for example, 
2.0E-6 represents a deadtime of 2 x 10-6 sec. If 
dead time is to be calculated from experimental 
data, TAU is set equal to zero. Dead time may be 
calculated for, say. Element 3 and then this 
value of TAU used also for Element 5. In this 
case, TAU for Element 5 would be set equal to 
3.0. The computer knows that this is not an ac- 
tual TAU value, because it expects the actual 
values to be less than one. A negative value for 
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TAU would indicate that no deadtime correction 
is to be made for that element. (We suggest that 
a probe analyst who chooses to regard deadtime as 
negligible is being overly optimistic.) If through 
some mischance the calculated value of TAU turns 
out to be negative, no deadtime correction will be 
made. 

2. Drift Correction 

Accurate microprobe analyses cannot be expect- 
ed unless the Instrument is operating in a stable 
manner with a minimum of drift. Program EMX cor- 
rects for drift by assuming a linear change with 
time between successive readings on the same 
standard. This correction is made with no regard 
to the amount of drift, and it is left to the dis- 
cretion of the operator to decide whether the 
drift rate is low enough to justify correcting and 
using the data. 

Successive data points on a standard were used 
to monitor drift during probe operation. Between 
each of these data points for the individual spec- 
trometers .Program EMX corrects the data for the 
variation in count rate on the standard by assum- 
ing drift to be a linear function of time. The 
elapsed time is measured by the accumulated count- 
ing time, and the time between successive readings 
is assumed to be negligible. If the elapsed time 
between successive readings is large, a Break Card 
is inserted into the data at this point. Drift 
corrections will then be made independently for 
the data before and after the break. No drift 
correction is made in the case where the variation 
between successive data points on the standard is 
within the deviation found in the successive points. 
In such a case Program EMX will print out a message 
indicating that no drift correction has been made. 

3. Background 

Program EMX incorporates the correction for 
background by two alternate methods. The method 
most suitable for the particular problem is select- 
ed by the probe operator, and this choice indicated 
by the appropriate symbol for MBKGD in the lead 
cards to Program EMX. Background may be measured 
in a single set of readings for the standard and 
another set for the sample, and an one scheme 
(MBKGD = 1) these values are subtracted from each 
standard and sample reading for the remainder of 
the problem. Further background (kt a points may 
be taken at any stage in the probe operation, in 
which case the new background values are adopted 
and subtracted from succeeding peak values. 

For precise analytical work and particularly 
for the measurement of minor elements in the range 
50 to 10,000 ppm it is necessary to monitor the 
background more carefully than can be achieved by 
the above method. In these cases background is 
measured for each data point on both sample and 
standard. In the second background method (MBKGD 
= 0) each peak value must be accompanied by a pair 
of background data points. The program averages 
these and subtracts the average value from the 
peak value. The sequence of background and peak 
values need not be kept constant (for example, see 
Fig. 2) , but two background values must be assoc- 
iated with each peak value. If one background 
reading is omitted or an extra value included, the 
wrong pairs of values will be averaged. 


Background values can be measured at the same 
wavelength as the peak values on a material simi- 
lar in composition to the sample (or standard) but 
lacking the element being analyzed. In the absence 
of such material it is convenient to measure back- 
ground data points at wavelengths above and below 
the peak wavelength. In practice, a series of data 
points are taken over a range of wavelengths in the 
neighborhood of the peak wavelength. The values 
of counts versus wavelength are examined to deter- 
mine the variation in peak and background intensi- 
ties and to search for any other emission peaks 
in this range. On the basis of these data, two 
wavelengths on either side of the peak are chosen 
(X + A, X - A) and pairs of background measure- 
ments made at these wavelengths. In the EMX pro- 
gram these pairs of data points are averaged to 
give the background values used. 


C. Output 

The final printed output from Program EMX 
consists of a tabulation of the concentrations 
of elements in each sample or sample point plus 
a tabulation of the elements measured, combined 
with appropriate s.d. components , as oxides, for 
example (see footnote page 4) . The element con- 
centrations, of course, are the first approxima- 
tions; i.e. , intensity ratios times the concen- 
tration of the element in its standard. Addition- 
al printed output includes all input data; input 
data after readings at each data point were aver- 
aged, after deadtime correction, and after cor- 
rection for drift; and the final measured x-ray 
intensities in counts per second with statistical 
errors. Figs. 5-9 are examples of printed out- 
put. Some time and paper can be saved by elimin- 
ating the intermediate output, but it is often use- 
ful to observe the effects of each correction as 
it is applied. 

At the option of the program user, the compu- 
ter will plot a graph of the concentrations (see 
Fig. 16). The specific plotting routines used in 
preparing this plot are usable only with a Calcomp 
165 incremental plotter used in conjunction with 
CDC 3600 and CDC 160-A computers; however, it 
should be relatively easy for laboratories with 
other types of computer-plotter facilities to 
substitute their own plotting routines into the 
program. 

When both programs are run together as a mul- 
tiple execution job, element concentrations from 
EMX are written on a scratch tape for input to 
Program EMX2. If data is to be used for EMX2 
input at some other time, it may be recorded on 
magnetic tape or punched cards. 
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III. PROGRAM EMX2 


A. Basic Correction Procedures 


a> 


I (A) b A 


( 6 ) 


Program EMX2 uses an iterative procedure to 
calculate composition from the input concentrations, 
which are first approximations calculated from mea- 
sured x-ray intensities by Program EMX. The neces- 
sary absorption and fluorescence corrections are 
made on the basis of the original input concentra- 
tions and a new composition obtained. This proce- 
dure is repeated as many times as specified by the 
program user or until the percentage change in con- 
centrations from one iteration to the next is less 
than a specified amount. 

The absorption correction used Philibert’s 
(3) formula for f(x)s 


f( X ) = 


l_+h 

u*x /cr )[ 1+ti 0+x /cr )] 


(3) 


where h = 1.2 A/=2, A and Z being the mean atomic 
number, x = u/p csc9, u/ p being the mass absorp- 
tion coefficient of the target for the measured 


v - — *b— a- , the effective o defined by 

E o • “ E c Duncumb and Shields (4) , E 

being applied voltage in kv and Ec the critical 
excitation potential for the particular element 
m kv. 


The fluorescence correction formula, a com- 
bination of those of Castaing (2) and Wittry (1), 
is as follows (in the case of Element A being 
fluoresced by Element B) : 


1(A) is the measured intensity of element A in its 
standard, I A is the measured intensity of element 
A in the sample, and S A is the concentration of 
the element in its standard. The correction factor 
R, which depends on the concentrations of all N ele- 
ments in the sample, is given by 

l £ 

R(C lt ...,C N ) = F (x)(l + ^) (7) 

The notation Rg ' refers to the correction factor 
for the standard used for Element A. 

In Program EMX2 BjuSls solved for the concen- 
trations Cj by the Gauss -Seidel iteration method, 
which has produced rapid convergence for all cases 
encountered so. far. The second-approximation con- 
centrations Cj 'are calculated as follows: The 

concentration ^ 0 f element 1 is calculated us- 
ing a correction factor R evaluated with the first- 
approximation concentrations CjW. The concentra- 
tion of element 2 is then calculated using a 

correction factor R evaluated with first-approxi- 
mation concentrations for elements 2 through N, 
but with the second-approximation value for Ci . The 
procedure continues in this manner, with R always 
using the best available (i.e. the most recently 
calculated) values for the C's. For the second 
approximation this procedure is expressed by the 
formula 


( 8 ) 


If A _ °>x B r / r-1 \ A /*b U B ln U B -U B -H 

I A 2 UB ‘V r /* B ^ompie U fl lnU A -U A + l 


In(H-u) + In(H-v) 


(4) 


where I is the emergent primary intensity of the 
characteristic wavelength of element A from direct 
excitation, 1^ is the emergent secondary intensity 
of the same wavelength from fluorescence by ele- 
ment B, W® is the k- fluorescence yield, Cg is the 
concentration of element B, r is the absorption 
jump ratio, A and B are the atomic weights of ele- 
ments A and B, yjj (^sample) mass absorption coef- 
ficient of the characteristic radiation of B in 
A (in the sample) , E^ and E® are the critical ex- 
citation potentials of elements A and B, - 
E q /eA, Ur = E 0 /E®, u = y |ample /tJ |ample sin e = 
XA/b| amp ^ e . and v = o/p| ain P le . 


The basic correction formula incorporating ab- 
sorption and fluorescence corrections is 


r - c (1) R S<A). 
C A " C A R 


(5) 


where is the first approximation concentra- 

tion of element A, 


This method would be prohibitively time consuming 
for hand calculation, but is natural for computer 
calculation. 


B . Physical Constants 


Data necessary for the calculation of the var- 
ious constants pertaining to absorption and fluor- 
escence corrections involving almost any elements 
are stored by Subroutine CONSTANT. Below are list- 
ed the data, sources of data, and methods of calcu- 
lating the required constants: 

1. SIN THETA. Sine of the takeoff angle, which 
with our microprobe is 52.5°. 

2. Atomic weights (A). From reports prior to 
1960 of the Commission of Atomic Weights, Inter- 
national Union of Pure and Applied Chemistry. 

3. Characteristic wavelengths (LAMBDA). For 
each element we have included only one character- 
istic wavelength, the one which we customarily 
use for analyzing the particular element. For 
elements with atomic numbers 5-35, LAMBDA is the 
weighted Ka wavelength (Kc^ given twice the weight 
of K 02 ); for atomic numbers 36-82, Laj; and for 
atomic numbers 83 to 90-92, MS. Wavelength data 
are from Bearden (7). 



4. Mass absorption coefficients (MU) are cal- 
culated by Heinrich’s (5) method: MU = CA n where \ 
is the wavelength of the radiation being observed 
and the coefficients (C) and exponents (EXP) are 
those determined by Heinrich. Values for absorp- 
tion edges (EDGE) are also from Heinrich's tables 
(5) . It should be noted that the data are not 
complete for all elements. As of now. Program 
EMX2 cannot be used for the analysis of extremely 
light elements (atomic numbers less than 12 ) . 

Values for C and EXP are also unavailable for the 
regions between the and absorption edges 
and beyond the N^ edge. The computer will calcu- 
late mass absorption coefficients as if the absorp- 
tion curves for these regions were merely continua- 
tions of those for the M 3 to M 4 and M 5 to regions. 
This is a very rough approximation and will intro- 
duce greater than usual errors when analyses in- 
volve elements such that the characteristic wave- 
length of one falls within the unknown absorption 
region of the other. 

5. Absorption jump ratios (R) used in the 
fluorescence corrections are calculated from 
Heinrich's data (5) by the formula R = (Ck/Cj.^) 

Ak ( n k " n kl^ > where C and n refer to the absorp- 
tion curve of the element fluoresced and = the 
characteristic wavelength of the fluorescing ele- 
ment. The absorption jump ratio is commonly cal- 
culated at the absorption edge, but this unneces- 
sary approximation is not made by Program EMX2 
(see Compton and Allison( 8 ), p. 542, or Philibert 
(3), p. 55). 

6 . K- fluorescence yields (OMEGA) are calcula- 
ted from the formula 

W 

( ■ ; ■ k 7 T ) 1/4 = A + BS + ca 3 (9) 

1 - w k 

using the coefficients of Hagedoorn and Wapstra(9), 
A =-0.064, E « 0.034, C = -1.03 x 10^>._ Because 
acceptable values of the L-fluorescence yields and 
ratios of K to L intensities are not easily evalu- 
ated at the present time (Philibert(3))*, only K-K 
fluorescence is considered. Corrections for K-L, 
L-K, and L-L fluorescence can be added to the pro- 
gram if reliable data become available. 


Except for the few introductory data cards (see 
below and Appendix 2) , the input to Program EMX2 
consists of the punched card or magnetic tape out- 
put from Program EMX. If no more than nine ele- 
ments are involved in each problem, there will be 
one card (or tape record) for each sample or point 
on the sample, and Program EMX2 can be run with 
Program EMX as a multiple-execution job. For 
samples involving 10 to 18 different elements, 
however, the data must be run through the first 
program as two separate problems, one including 
nine elements and the other including the remainder. 
The punched card output must then be put together so 
that the two cards for each sample are together in 
the data deck for EMX2. 

Sometimes samples contain elements such as car- 
bon that cannot be analyzed directly by many micro- 

*There is some difference of opinion about this; 
see Reed (10) . 


probe facilities. Concentrations of such elements 
must be considered in the total composition, how- 
ever, if absorption and fluorescence corrections 
are to be made. The probe operator can estimate 
the concentrations of elements which were not or 
could not be analyzed and add them to the data 
cards for EMX2 after the concentrations of the 
other elements, setting NFIX on the Identification 
Card (see Appendix 2) equal to the number of such 
elements. Concentrations cf these fixed elements 
will be used in calculating average atomic number, 
etc. , but will not themselves be corrected. 

For a particular element i the correction factor 
R(i) combines the absorption and fluorescence cor- 
rections for the standard (see section I II A) . The 
value of R(i) depends on the composition of the 
standard and the operating voltage. To solve 
Equation 5, Program EMX2 will calculate R* for 
the unknown, but the values of R(i) , the correc- 
tion factors (SCF) for the elements in their stan- 
dards, must be read into the program on the Ele- 
ment Cards (Appendix 2) , Thus before any analyses 
of actual samples are corrected, the necessary 
correction factors must be calculated. This is 
done by considering the standard as a separate 
problem, with an Element Card for each element of 
the standard on which SCF is set equal to zero. 

The operating voltage (EO) for this problem is 
set equal to that actually used in the analyses. 

The computer will calculate each SCF value and 
punch new Element Cards containing all the nec- 
essary information including the SCF. These 
cards may be saved for use as Element Cards in 
any problem in which the same standard and E0 
were used. The most efficient procedure is to 
build up a file of Element Cards for all the 
standards commonly used (pure element standards 
as well as mineral or alloy standards) with cor- 
rection factors evaluated for all E0 which are 
likely to be used. 


D. Computer Method 

After the preliminary calculations of the 
physical constants (listed under Sec. IIIB, above) 
by Subroutine SETUP, the computer looks at the 
first SCF (see Sec. IIIC) to see if the correc- 
tion factors for the standard are known. If they 
are not (i.e., if SCF = 0), the problem is assumed 
to be one of determining the correction factors 
for a particular standard. In this case the pro- 
gram calculates the SCF, punches new element cards, 
and proceeds to the next problem. If the SCF do 
not equal zero, the computer proceeds to the sample 
data and reads in the first approximation concen- 
trations (OCON(J), J = 1 to NEL, where NEL = number 
of elements analyzed) for one sample or point on 
the sample. Subroutine OXYGEN is called to calcu- 
late the concentrations of oxygen (OCON(NELPI) ) and 
silicon (OCON (NEPL2)) (if 0 or SI factors present) 
as a function of the concentrations of the other 
elements. Concentrations of all elements and their 
total are printed, and the computer then proceeds 
to the corrections themselves. 

The first step is the storage of all OCON m 
CON and the multiplication of each OCON by its SCF. 
Now OCON is not the first approximation but a pro- 
duct of the first approximation and the R correc- 
tion factor for the standard according to Equation 
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5. This product will remain constant throughout 
the iterations. 

FCHI, the absorption correction factor, is 
calculated in Function ABSORB. 

Function FLUOR(NN) calculates the fluorescence 
correction term, which may include the fluorescence 
effects of from zero to four other elements. Sub- 
routine SETUP has previously determined which ele- 
ments fluoresce each of the elements analyzed in 
the problem. For simplicity, the former may in- 
clude only the four elements nearest in atomic 
number whose Ka lines are less than the K absorp- 
tion edge of the element in question. Thus, for 
example, in determining which elements might fluor- 
esce Fe (Z - 26) , the program looks only for ele- 
ments in the problem with atomic numbers in the 
range 28-31. Although fluorescence by elements 
with higher atomic numbers or by the KfJ line of 
Co (Z = 27) might be significant, these are not 
taken into consideration. 

In the Fortran notation of the computer pro- 
gram, the function calculated by Function FLUOR (NS) 
is 

NF(NN) 

FLUOR - £ COU L '.Vj - 

3 • 1 


where L is index of fluorescing element, NN is in- 
dex of element fluoresced, NF is number of elements 
fluorescing NN, FMUBAR is mass absorption coeffi- 
cient of alloy for wavelength of element L, MUBAR 
is mass absorption coefficient of alloy for wave- 
length of element NN, 

mr = In (1 + MUBARAFMUBAR-S INTHETA) 
MUBAR/(FMUBAR • S INTHETA) 


In (1 + SIGMA/ FMUBAR) 

S I GMA/ FMUBAR 

and U is the U from Subroutine SETUP, a variable 
which contains all the parameters of the fluores- 
cence correction equation except those that de- 
pend on the concentrations of elements. Note: U 
here is not just EO/EC. 

A new approximation (CON 2) is then calculated 
according to Eqn. 8, which in Fortran notation is 
written as C0N2 = OCON(I)/ (FCHI* (1.0 + FL)). 

DELCON , the change in the concentration of Element 
1 that resulted from this formula, is calculated; 
concentrations of s.d. components (if any) are 
modified, and the value of C0N2 replaces CON(I). 
Then the procedure is repeated for all the other 
elements. Note that the estimated composition, 
from which FCHI and FL are calculated, changes 
with the correction of each element; the same 
estimated composition does not carry through an 
entire iteration. The results of each correction 
are printed out (see Fig. 13), and the total of 
concentrations is printed after each iteration. 

The limits of the iteration have been read 
into the computer from the Problem Card, with 
IMIN and IMAX specifying the minimum and maximum 


number of iterations to be performed. Within 
these limits, the correction procedure will be 
iterated until DELCON for any element present in 
amounts greater than one percent does not change 
by more than DMIN from one iteration to the next. 

No method for calculating estimated errors in 
the final results is included in the program. 


E. OUTPUT 

Final output for Program EMX2 is in the same 
format as that from Program EMX, including the 
optional graph. (If more than ten elements are 
included in the problem, there will be two graphs 
because too many lines on one plot would be un- 
decipherable.) 

Concentrations in weight percent shown for 
each element and the total for all the elements 
analyzed are those calculated by the correction 
procedures described above. They are not normal- 
ized to sum to 100 percent. The information that 
the total corrected concentrations for a particu- 
lar problem vary quite markedly from 100 percent 
can give a useful indication of how well the micro- 
probe was working, how carefully the analyses were 
performed, or perhaps how the computer correction 
method needs to be revised. We do not believe 
that this valuable information should be lost. 

Printed output from EMX2 begins with the input 
data from the Identification and Element Cards as 
well as all the mass absorption coefficients, wave- 
lengths, etc., determined by Subroutine SETUP. 

Next comes one printed page for each sample or 
point on the sample showing the first approxima- 
tions (those read in), and the concentrations, 
f(x), and fluorescence correction values calcu- 
lated in each iteration, plus the change in the 
concentration of each element from one iteration 
to the next. Figs. 11 to 16 are examples of 
printed output from Program EMX2. 
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Standard for Elements 1, 2 and 
3 on spectrometers 1, 2 and 3, 
respectively. NAME indicates 
that the mineral standard used 
was Marjalahti. Identifying 
names are optional; this col- 
umn could have been left blank. 

Background readings for stan- 
dard } above and below peak. 

Sample background readings. 


Sample readings on all three 
spectrometers. 


Standard readings repeated. 
Computer will interpolate 
between this and previous 
group of standard readings 
to correct for drift. 

"BREAK" Card: indicates end 
of a data set. 

Standard peak for elements 5 
and 6 on spectrometers 2 and 
3, respectively. Zero in the 
second column of the KODE 
indicates that reading on 
spectrometer 1 is to be 
ignored. 
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Standard for Element 4 on 
spectrometer 1; other spec- 
trometers ignored. 

First sample background read- 
ing. Will be averaged with 
the one that comes after the 
sample readings. Note that 
pairs of background readings 
need not appear consecutively, 


999999 Card signalling the end of the 

data for this problem. 


Figure 1. Data sheet as it comes from the microprobe illustrating the use of the six-digit code. The 
columns, from left to right, are counting time in hundredths of a second, sample current in pa x 10“4, 
counts on spectrometers 1, 2 and 3, the six-digit code (KODE), and, if desired, a name for identification 
(NAME). In this problem NS = 1, so each reading on the sample is a separate data point. MBKGD = 1 so 
background data points do not accompany each individual peak data point. 
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000100 000800 003837 000000 000000 410000 Deadtime data on Spectrometer 

001000 000780 003712 000000 000000 410000 1. Printout columns for spec- 

001000 000760 003654 000000 000000 410000 trometers 2 and 3 were not 

used, so turned off. NDT - 1, 
so each reading is one data 
point. Of course many more 
dead time data points would 
be taken. 


"BREAK" Card. Numbers which 


000000 

000013 

000476 

000000 
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000000 

happened to get printed in 
cols. 2 and 3 will be ignored. 

002000 

000000 
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000000 
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Standard peak and background 
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that readings in this problem 

010000 
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110000 

were taken over different time 
intervals. 
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000000 

110000 

First sample. MBKGD equals 0, 
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000000 
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110000 

so a pair of background data 

010000 

000000 

000695 

000000 

000000 

010000 

points will be taken with each 
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Second sample. Note that back- 
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ground points were taken before 
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the first sample but afte** this 
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sample. Computer will keep 
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track of which background goes 
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with which sample. 
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A standard peak data point taken 
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000000 
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000000 

000000 

210000 

again to check for drift. Be- 

002000 
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000000 

000000 

210000 

cause MBKGD equals 0 but no 

002000 
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118216 

000000 

000000 

210000 

background readings were taken 
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000000 
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110000 

for this standard peak, this 
data point will not be cor- 

010000 

000000 
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000000 

000000 

110000 

rected for background. But 
that doesn't matter because 

010000 

000000 

000708 

000000 

000000 

010000 

010000 

000000 

000721 

000000 

000000 

010000 

this point is used only for 

010000 

000000 

000501 

000000 

000000 

110000 

the drift correction. 

010000 

000000 

000519 

000000 

000000 

110000 


000000 

000000 

000000 

000000 

000000 

999999 



Figure 2. Data as it comes from the micrpprobe illustrating the use of the six-digit code. First three 
readings are part of a data set used to measure dead time. The remaining readings refer to a problem in 
which NS = 2 so that each standard background, sample peak and sample background data point comprises two 
readings which will be averaged. MBKGD =0, so each peak data point is accompanied by a pair of background 
data points. 
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9 .30436205-0? 


IQ . 

1. 63540005 0 0 

5.7161§565 

.0? 

5 S 6978286E 

02 

-1 . 83570 0 3E 

no 

1.14439955-01 


11 

2.06980005 on 

5,72893525 

n? 

5 1 . 6910365E 

0? 

**3 . 7898709E 

00 

1,4451 55 IE- 01 


12 

1.624 3 nnoF on 

5.676?07QF 

n? 

5 1, 69800225 

0? 

2.1095222F 

no 

1.14449655-01 


13 

l . 968 7 n n o E no 

5.61 00472E 

02 

5 , 6926172E 

02 

8.2569992F 

00 

1 ,40369565-01 


14 

2 t 51 755 0 0E no 

5,66999665 

D? 

5. 684Hi555F 

„Q2 

1 . 40389O5E 

00 

1.77605045-01 


15 

3.3996500F no 

5.691 7624F 

02 

5.6702432F 

0? 

-2.1519188F 

nn 

2.3891 721F-01 


FOR ELEMENT 1? ON SPECTROMETER 3* TAlj’ = ' 4.775-06 SECONDS 


Figure 4, Output of dead time calculation from 
Program EMX. 



PROBLEM NO, 


IB 


S^ESSARTINE RUTHBRFORO 3 EGMATITE ARCH MAR 1966 


NO, OF ELEMENTS = 5 
0 DEADTIME READINGS PER S = r 


2 SAMPLE 

READINGS 

PER SET 





background reading 

method = 1 





SIGMA FOR 

STANDARD 

peak mat not 

■ EXCEED 3, 

,0 * 5QR 

TF<N) 



ATOMIC PER CENT A 

TAJ 

oxygen 

SILICON 


element 

number of 

■ standard 


factor 

factor 


1 

25 

34, 970000 

1 .60E-05 

,29110 

- 0 

MNO 

2 

20 

6,980000 

2, 00F-06 

,39220 

-o 

CAO 

3 

13 

14 . 700030 

2 , OOE-06 

. 88960 

- n 

A l_ 2 O 3 

4 

26 

6 . 70 0 0 0 0 

1 , 60E-05 

, 28650 

-o 

FEO 

5 

14 

IB. B0 00 00 

2, 00E-06 

1,13900 

-0 

S I 02 


Figure 5. Output from Program EMX. First output 
page showing problem description. 
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PROBLEM NO. 18 S=>lr 5S ART t NE RUTHERFORD PEGMATITE ARCH MAR 1966 


IvjPUT data 


NAME 

TIME CURRENT 

SPEC 1 

SPEC 2 

SPEC 3 

KODF 

1 

20.00 

0 

82609 

27242 

53 

210000 

1 

20.00 

0 

83333 

27554 

59 

21 0000 

1 

20 . 00 

0 

81919 

27580 

54 

210000 

1 

20,0 0 

0 

82353 

27778 

62 

210000 

1 

20 . no 

0 

82898 

27674 

67 

21 0000 

2 

20,00 

0 

29*5 

60848 

18813 

202223 

2 

20,00 

0 

280 

60669 

18580 

202223 

2 

20.00 

0 

237 

60744 

18683 

2022?3 

2 

2 0,00 

0 

313 

60603 

18996 

2022?3 

2 

20 .00 

0 

287 

59902 

185 94 

2 0 22?3 

3 

20, 00 

n 

45635 

11379 

11872 

010 203 

3 

20,00 

n 

45679 

11379 

11649 

Q 1 0 2 0 3 

4 

20.00 

n 

45231 

11307 

11789 

010203 

4 

20.00 

0 

4 5221 

11314 

118? 0 

0 1 0 2 0 3 

5 

20 .00 

0 

4 4 8 91. 

10515 

11886 

010203 

5 

20.00 

n 

45595 

10525 

11974 

010 2 0 3 

6 

20,00 

n 

44192 

12028 

11376 

0 1 020 3 

6 

20.00 

n 

44215 

12215 

11977 

0 1 0 2 0 3 

7 

20.00 

n 

450 40 

12840 

11903 

0 1 0 2 0 3 

7 

20.00 

0 

450 22 

12616 

11756 

010203 

8 

20.00 

0 

39273 

14946 

U8?5 

01 0?03 

8 

20.00 

0 

39605 

15069 

1 18?B 

0 1 0 20 3 

9 

20,00 

0 

44900 

12495 

11831 

01 0 2 0 3 

9 

20.00 

n 

44371 

12478 

11799 

010203 

10 

20.00 

0 

4 63 59 

11800 

12004 

0 10 2 0 3 

10 

20.00 

n 

46724 

11691 

11951 

010203 

11 

20.00 

n 

45120 

12941 

11895 

0 1 n ? 0 3 

11 

20.00 

0 

450 04 

12911 

11780 

010203 

12 

20.00 

n — 

“52604 

10574 

11695 

01 0203 

12 

20.00 

0 

52448 

1 0457 

11767 

01 0203 

13 

20.00 

0 

270 

61696 

18608 

202223 

13 

20.00 

0 

302 

60964 

18670 

202223 

14 

20.00 

0 

81351 

27847 

76 

21 2020 

14 

20.00 

0 

80964 

? 7808 

51 

212020 

15 

20,00 

0 

452 7? 

11313 

119 3 0 

010203 

15 

20.00 

0 

45801 

11459 

11604 

010203 

16 

20, GO 

0 

54102 

10271 

11781 

010203 

16 

20 .00 

0 

54381 

1 0 0 70 

11 9?6 

010203 

17 

20.00 

0 

520 73 

10510 

11930 

010203 

17 

20 .00 

0 

5150? 

10468 

11615 

010203 

18 

20 .00 

0 

52731 

10313 

11464 

010203 

10 

20 .00 

0 

52212 

10341 

118?6 

010203 

19 

20,00 

n 

52819 

8888 

11689 

010 20 3 

19 

20 .on 

n 

53108 

9048 

11855 

010203 

20 

20.00 

n 

52501 

10685 

11628 

010203 

Figure 6 . Output 
identification of 

from Program EMX. 
data points. 

Second 

output page 

showing 

raw probe data and code 



PROBLEM NO, 


18 


PAGE 3 


SafcSSARTINE RUTHERfORO 3 EGMATIT= ARCH MaR 1966 

a vruiEinrm ~ 

REL. SPECTROMETER 1 SPECTROMETER 2 SPECTROMETER 3 


1 

time 

50 

current 

9 

CTS/SEC 
4131.1 _ 

SIGMA 

10,7 

2 

1 

CTS/SEC 

0 

SIGMA 

0 

0 

0 

CTS/SEC 
0 ’ 

SIGMA 

0 

0 

0 

2 

150 

' o' 

n 

"0* 

2 

0 

30 277 7 

7,5 

2 

2 

936.7 

3.5 

2 

3 

3 

220 

0 

2284,1 

7,6 

0 

1 

569. 0 

3.0 

0 

2 

588,9 

3.9 

0 

3 

4 

260 

0 

2261.3 

7,5 

0 

1 

565 • 5 

3,8 

0 

2 

590,2 

3 . 8 

0 

3 

5 

300 

0 

2262. 2 

12.4 

9 

1 

526.0 

3.6 

0 

2 

596.5 

3.9 

0 

3 

6 

340 

a 

221(1 , 2 

7,4 

0 

l 

604,1 

3.9 

0 

? 

596,3 

3.9 

0 

3 

7 

380 

0 

2251.6 

7,5 

0 

1 

636. 4 

4,0 

0 

2 

591,5 

.3.8 

0 

3 

a 

420 

a 

1951.9 

0,3 

0 

1 

750.4 

4.3 

0 

2 

591.3 

3.0 

0 

3 

9 

460 

0 

2231.0 

9,4 

0 

1 

624.3 

4.0 

0 

2 

590.8 

3.8 

0 

3 

10 

500 

0 

23"27 , 1 

7,6 

0 

1 

507.3 

3.8 

0 

2 

598,9 

3.9 

0 

3 

11 

540 

3 

2253.1 

7,5 

0 

1 

646.3 

4,0 

0 

2 

591.9 

3.8 

0 

3 

12 

580 

n 

2626.3 

8,1 

9 

1 

525.8 

3,6 

0 

2 

586.5 

3.8 

0 

3 

13 

620 

0 

fi 

0 

2 

1 

3064.5 

12.9 

2 

2 

932.2 

4.8 

2 

3 

14 

660 

0 

4050,4 

10,1 

2 

1 

0 

0 

2 

0 

n 

0 

2 

0 

15 

700 

3 

2276,8 

9.4 

9 

l 

569. 3 

3.8 

o 

2 

588.4 

5,8 

0 

3 

16 

740 

3 

2712,1 

8,2 

0 

1 

503.5 

3,6 

0 

2 

592.7 

3.0 

0 

3 

17 

700 

a 

2589 , 4 

10,1 

0 

1 

524.5 

3.6 

0 

? 

583,6 

5.6 

0 

3 

10 

820 

0 

2623.6 

9,2 

9 

1 

516 • 4 

3.6 

0 

2 

582,3 

6 . 4 

o 

3 

19 

3h0 

0 

2640.? 

8,1 

a 

1 

4 4A.4 

3.3 

0 

2 

588.6 

3,8 

0 

3 

20 

903 

0 

2622.5 

8,1 

0 

i 

53?. 3 

3.6 

0 

2 

578.9 

3.8 

0 

3 

21 

940 

o 

2735.7 

6,3 

0 

l 

516.7 

3.6 

0 

2 

581.3 

3.8 

0 

3 

22 

980 

3 

2670.5 

8.2 

0 

l 

519 • 4 

5 . 4 

0 

2 

581.0 

3.8 

0 

3 

23 

1020 

0 

2650 , 7 

IP ,0 

0 

l 

472-5 

3.4 

0 

2 

574.1 

3.3 

0 

3 

24 

1060 

0 

2995.2 

14.6 

9 

l 

417.4 

3.5 

0 

2 

581.7 

3.8 

a 

3 

25 

1100 

0 

n 

0 

2 

i 

2950.7 

3.6 

2 

2 

899,6 

4.7 

2 

3 

26 

1140 

0 

3939.9 

1C, 4 

2 

l 

0 

0 

2 

0 

o 

0 

2 

0 

27 

1180 

0 

37.8 

1,0 

3 

i 

0 ’ 

0 

0 

0 

9 

0 

0 

0 

28 

122*1 

0 

32.5 

. 9 

3 

t 

0 

0 

0 

0 

9 

0 

0 

0 

29 

1260 

0 

r 

0 

3 

0 

64 . 7 

1.3 

3 

2 

12.1 

.5 

3 

3 

30 

1300 

0 

0 

0 

3 

0 

64.5 

1.3 

3 

2 

18,3 

, 7 

3 

3 

31 

1340 

0 

25. r 

,8 

1 

1 

70.6 

1.4 

l 

2 

11.8 

.5 

1 

3 

32 

1380 

0 

25,9 

1,3 

1 

t 

64.7 

1.3 

1 

2 

9.3 

.5 

1 

3 

33 

20 

0 

4156.3 

10,2 

2 

l 

0 

0 

2 

0 

0 

0 

2 

a 

34 

60 

0 

n 

0 

2 

3 

3101.6 

8.8 

2 

2 

945.2 

5,0 

2 

3 

35 

103 

n 

3066.5 

8.8 

n 

1 

447,4 

5.9 

0 

2 

608.9 

3.9 

0 

3 

36 

140 

3 

3212.4 

10,8 

0 

1 

413,1 

3.2 

0 

2 

598.5 

.3.9 

a 

3 

37 

18 n 

0 

3-121.5 

8,7 

0 

1 

467.3 

4.8 

0 

2 

60 6.5 

3.9 

0 

3 

38 

220 

0 

3139.2 

6,9 

0 

l 

42?, 7 

3,3 

0 

2 

599.9 

3.9 

0 

3 

39 

260 

n 

?°72 . 9 

8.6 

8 

1 

464.5 

3.4 

0, 

2 

599, 9 

3.9 

a 

3 

40 

30 0 

0 

3306.5 

21,8 

0 

t 

374.0 

3.1 

0 

2 

605.7 

3.9 

0 

3 

41 

340 

3 

3342.1 

9,1 

9 

1 

384.4 

3.1 

0 

2 

603.0 

3.9 

a 

3 

42 

3du 

0 

3?97.? 

9.0 

9 

1 

394.6 

3.1 

0 

2 

596.2 

3.9 

0 

3 

43 

420 

0 

3215.5 

9.0 

0 

1 

333.1 

3,1 

0 

2 

594,9 

3.9 

0 

3 

44 

460 

0 

3249.6 

9,0 

0 

1 

349.1 

3.0 

0 

2 

583.9 

3.8 

0 

3 

45 

50 u 

0 

9 

0 

2 

3 

3033*1 

8.7 

2 

2 

936.1 

4 . 8 

2 

3 

46 

540 

9 

4134.7 

10.2 

2 

1 

0 

0 

2 

0 

0 

0 

2 

0 

47 

' 500 

0 

3220.5 

19,4 

i) 

1 

340.3 

3.0 

0 

2 

591.9 

3.8 

0 

3 


Figure 7. Output from Program EMX. Averaged data. Columns from left to right are data point number, rela- 
tive time (for drift correction), sample current (not measured in this problem), averaged count rate in t-> 

counts per second for each data point on spectrometer 1, error, code number showing type of reading, element 
number, averaged count rate on spectrometer 2, etc. 
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PROBLEM NO, 18 S=>£3SARTINE RUTHFRF3R0 PEGMATITE ARCH MAR 1966 

CONCENTRATIONS CORRECTED FOR D5ADT I ME j DRIFT, AND BACK 

SAM?lfc CORRECTED CQNCEN* RELATIVE GROUND 

E sample current CTS/SEC tration ERROR 


1 1 STD 0 

3 1 0 

4 1 0 

5 1 0 

6 1 0 

7 1 0 

8 1 0 

9 1 0 

10 1 0 

11 1 b 

12 1 0 

15 1 0 

16 1 U 

17 i 0 

18 1 0 

19 1 0 

20 1 0 

21 1 0 

22 1 0 

23 1 0 

24 1 0 


2 2 STD 0 

3 2 0 

4 2 0 

5 2 0 

6 2 0 

7 2 0 

8 2 0 

9 2 0 

10 2 0 

11 2 0 

12 2 0 

15 2 0 

16 2 0 

17 2 0 

18 2 0 


19 2 0 

20 2 0 

21 2 0 

22 2 0 

23 2 0 

24 2 0 


4388,3 

.34970 

,002 

2357.8 

, 13789 

t 004 

2336,1 

,19616 

, 0 0 4- 

2339 , 9 

.19647 

, 0Q6 

2286,5 

.19221 

,004 

2334 , 3 

.13602 

,004 

2012,7 

, 16039 

,005 

2318,7 

.13477 

,005 

2425,7 

, 19330 

,004 

2347,8 

,18710 

,004 

2761.7 

.22007 

,004 

2343,6 

.13676 

,005 

2824.5 

,22508 

,004 

2697. 1 

.21493 

,004 

2741 . 8 

,21849 

.004 

2776.3 

,22124 

,004 

2755 , 2 

,21956 

,004 

2888,5 

, 23018 

,004 

2823 , 5 

,22500 

,004 

2808 , 9 

, 22384 

,004 

3204 , 4 

. 25535 

,005 

2979,5 

, 0 6980 

,003 

498,3 

.01167 

,008 

494 , 3 

,01-158 

, 008 

454 , 2 

, 01064 

, 009 

533,5 

,01250 

,008 

563 , 1 

,01319 

,008 

675,8 

. 01583 

,007 

549,7 

. 01288 

,008 

512,2 

,01200 

,008 

570,1 

, 01336 

,008 

450 . 0 

. 01054 

,009 

499.8 

,01171 

,008 

438,9 

, 01028 

,009 

454 , 8 

, 01066 

,00b 

446, 7 

,01047 

,009 

378,6 

, 0 0887 

,010 

462 , 7 

, 01084 

,009 

447 . 1 

,01047 

,009 

440 . 8 

,01033 

,013 

402.7 

.00943 

,009 

347.5 

, 00814 

,011 


Figure 8. Output from Program EMX. First approximation. Columns from left to 
right are data point number, element number, sample current, count rate in counts 
per second corrected for deadtime, drift and background, corresponding composi- 
tion (first approximation) and statistical error. 
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PROBLEM NO. IB S»fc5SARTINE RUTHERFORD p tQMAT J TE ARCH MAR 1966 


CO\lCENTRATlONS CORRECTED FOR D5ADT I ME , DRIFT, 


total 

25 

20 

55, 1833 

18,7391 

1.1674 

56.5262 

18.6158 

1,1580 

56.6706 

18.6466 

1.0641 

57.3713 

18.2205 

1.2498 

57.0652 

18.6017 

1.3190 

58.3307 

16,0391 

1.5832 

57.8720 

18,4771 

1.2877 

58.2862 

19,3300 

1,1999 

58,2711 

18,7095 

1,3356 

57.0310 

22,0074 

1.0543 

57.2659 

18,6760 

1.1709 

58.0152 

22.5078 

1.0281 

57.6675 

21.4926 

1,0655 

57.3793 

21.8494 

1.0465 

57,0401 

22,1240 

. 8870 

57.9029 

21,9561 

1.0840 

58.2234 

23,0179 

1.0474 

57.4976 

22,4998 

1. 0326 

57.2470 

22.3340 

.9435 

57,8181 

25.5354 

.8142 

58.4975 

25.3499 

.8650 

59.0397 

26,6425 

. 7999 

59. 0839 

24.9774 

.9150 

58.7357 

26,0216 

. 8144 

58.7214 

24,5735 

.9129 

59,0574 

27.5203 

. 7055 

59,0784 

27.8477 

.7315 

58,1192 

26.6658 

, 7570 

57,7727 

26.7513 

.7319 

56,9297 

27.0646 

, 6544 

57,0592 

26.7337 

. 6362 

56.8160 

27,8768 

.5787 

56,4300 

26.0319 

, 6555 

56,2167 

27,3320 

. 5705 

56,2717 

26 .7507 

, 5634 

56,4949 

27,5216 

. 4786 

56.2432 

28,7223 

. 4996 

57,1427 

27,2230 

. 6346 

57,0323 

26,8648 

.6313 

56,9173 

28,5053 

,4643 


13 

26 

14 

9,2061 

8,8023 

17.2184 

9.2412 

10,1431 

17.3682 

9.3413 

10,3872 

17.2315 

9.3385 

11,2715 

17.2909 

9.2619 

10.6624 

17.2201 

9.2587 

14,0259 

17.4238 

9,2495 

11,4579 

17.3998 

9,3792 

11,163? 

17.2139 

9,2675 

11,5483 

17.4102 

9,1825 

7,4855 

17.3021 

9.2112 

10,6452 

17.5625 

9,2303 

7,6470 

17.5519 

9,2156 

8,, 3213 

17.5724 

9,1139 

7,7834 

17,5860 

9,2152 

7,3558 

17,4581 

9,0600 

8,2380 

17.5648 

9,0991 

7,4067 

17.6523 

9.0947 

7,4176 

17.4529 

8,9834 

7,3607 

17,5754 

9,1043 

4,6985 

17,6657 

9,4608 

5,1888 

17.6330 

9,3045 

4.5568 

17.7360 

9,4385 

6,0359 

17.7171 

9,3422 

5 , 0 ? 9 ? 

17,5282 

9.3518 

6.2747 

17.6085 

9.4517 

3,6528 

17,7271 

9 . 4180 

3 .4825 

17.5988 

9,3182 

3 , 9613 

17.4169 

P.3055 

3.7337 

1 7.2503 

9,1388 

2.8657 

1 7.2062 

9,1823 

3,0256 

17,4313 

9 , 0779 

2 , 0738 

17,2087 

9,2305 

3,1105 

17.4015 

9,0234 

2,0134 

17.2773 

8 , 8692 

2 , 8402 

17.2483 

9 , 0119 

1 , 9081 

17,5747 

8,9815 

, 8430 

1 7 . 1 969 

9,0060 

2.7017 

1 7,, 577 4 

8,9311 

3,2180 

17.3371 

9,0080 

1 ,2703 

1 7 , 6694 


AND BACK 
GROUND 


Figure 9. Output from Program EMX. First approximation, 
weight percent, elements identified by atomic number. 


Concentrations in 
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PROBLEM NO, 18 S=>fc 5SART I NE RUT HERE ORE) PEGMATITE ARCH MA& 1966 

CD\JCENTRAT 1 0N8 CORRECTED FOR DEADTIME, DRIFT, AND BACH 

GROUNT 


TOTAL 

MN0 

CAO 

AL 203 

FED 

S I 02 

91.4341 

24,2586 

1,6253 

17.3958 

11,3241 

36,8302 

93.3087 

24 ,0 348 

1,6121 

17 ,4621 

13 t 0490 

37,1505 

93 , 4286 

24 ,0 746 

1,4815 

17,6513 

13,3631 

36,8581 

94,3966 

23,5545 

1.740 0 

17,6461 

14.5008 

36,9853 

93,9054 

24,0167 

1.8364 

17,5013 

13,7171 

36,8339 

95.7213 

20 . 7081 

2,2041 

17,4953 

18 .0443 

37.2695 

95,0851 

23,8557 

1.7927 

17,4779 

14,7406 

37,21-82 

95.5324 

24,9570 

1.6705 

17,7230 

14,3614 

36,8205 

95,6244 

24.1559 

1,8594 

17.5119 

14,8568 

37.2404 

93,8720 

28,4137 

1,4677 

17,3513 

9,6301 

37,0092 

94.4096 

24,1126 

1 , 630 0 

17,4056 

13,6951 

37.5663 

95 .4088 

29, 0599 

1.4314 

17.5360 

9,8379 

37.5436 

94 , 9392 

27,7491 

1,4834 

17,4139 

10,7053 

37.5875 

94,5182 

26 ,2098 

1.4570 

17.2216 

10,0134 

37.6164 

94,0183 

26,5643 

1.2349 

17,4131 

9,4632 

37.3428 

95.1457 

28.3475 

1 .5092 

17,1198 

10,5981 

37,5710 

95.6573 

29,7184 

1,4582 

17,1937 

9,5288 

37,7582 

94 , 5470 

29, 0495 

1,4376 

17.1854 

9,5427 

37.331 7 

94.2519 

28.8999 

1,3135 

16,9751 

9, 4695 

37.5938 

95 .1373 

32,9638 

1,1335 

17,2035 

6,0446 

37,7869 

96,2030 

32, 7292 

1.2042 

17,6771 

6,6754 

3^.7170 

96,8932 

34,3931 

1.1136 

17,5318 

5,8623 

3 7 . 9374 

97 . 0193 

32,2484 

1.2738 

17.8350 

7.7651 

37,8970 

96.3463 

33,5965. 

1,1338 

17,6531 

6,4701 

37,4929 

96, 4060 

31 .7268 

1.270 9 

17.6712 

8, 0724 

37 , 6646 

96, 9912 

35.3314 

. 9822 

17.8599 

4,6993 

37.9184 

96,8-927 

35.9541 

1,0183 

17.7963 

4 , 4802 

37 . 6438 

95 . 4407 

34,4?82 

1,0539 

17,6077 

5,0962 

37.2547 

94,8430 

34,5386 

1.0190 

17.5836 

4,8033 

36,8984 

93 , 6136 

34,9430 

.9111 

17.2607 

3,6868 

36,8041 

93.9952 

34,5305 

,8858 

17.3509 

3,8925 

37.2856 

93.4285 

35 ,9918 

.8057 

17,1537 

2 , 6679 

36, 8094 

93,1880 

33,6098 

,9126 

17,4420 

4 , 0017 

37.2219 

92,6797 

35,2384 

.7943 

17,0506 

2.5903 

36.9562 

92,6294 

34,5378 

. 7843 

16,7592 

3,6539 

36.8942 

93,2753 

35,5331 

, 6664 

17.0289 

2,4548 

37.5922 

92,6188 

37,0933 

.6955 

16,9714 

1 . 0845 

36,7841 

94,1226 

35,1476 

.8835 

17.0177 

3,4758 

37,5981 

93.7587 

34.6351 

, 8789 

16,9706 

4.1400 

37.0840 

93. 9 0 0*3 

36,8032 

. 6465 

17,0214 

1,6342 

37, 7949 

94 . 6211 

34,2531 

. 9794 

17,7104 

4,6581 

37,0121 

94,7495 

37.0218 

, 7655 

17,6435 

1,4096 

37 .9090 

95,6673 

38,5120 

.6120 

17,8083 

,9866 

37.7484 

94,8211 

37,6485 

, 7903 

17,8053 

1,2288 

37.3481 

95,7229 

37,9347 

.7728 

17.9946 

1,2009 

37,8199 


Figure 10. Output from Program EMX. First approximation. Concentrations in 
weight percent, expressed as oxides. 
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PROBLEM NO, 18 S^SSARUNE RUTHERFORD PEGMATITE ARCH MAR 1966 

1, STANDARD FOR ELEMENT T)XYgEN SILICON 

NL Z CONCENTRATION F( CHI ) *( 1+FLUOR ) FACTOR FACTOR 


1 

25 

,349690 

2 

20 

, 069300 

3 

13 

,147000 

4 

26 

,087000 

5 

14 

,188000 


964476 

,291100 

-0 

MNO 

874575 

,399200 

-0 

CAO 

599711 

,889600 

* 0 

AL203 

974452 ’ 

,286500 

-0 

FEO 

528008 

1,139000 

*0 

S 1 02 


mass absorption coefficienTS---mu/rh'0 


EMITTER= 25 
ABSORBER 

20 

75 

79 

286 

20 

371 

139 

13 

117 

432 

26 

89 

321 

14 

145 

531 

8 

31 

116 


13 

26 

14 

3419 

64 

2225 

1669 

296 

1086 

■3B6 

93 

3493 

3844 

71 

2502 

"50 4 

115 

328 

1505 

25 

966 


Figure 11. Output from Program EMX2. First out- 
put page showing for each element an identifying 
number, atomic number, concentration of the ele- 
ment in its standard, the correction factor (R) 
for the 20 kv operating voltage (see Section IIIA) , 
oxygen factor, silicon factor, and formula of the 
oxide. Calculated mass absorption coefficients 
are shown in the lower table. 


20 K' 
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PROBLEM NO. 18 So=SSARTlNE RUTHERFORD PEQMAT I Tg ARCH MAR 1966 20 KV- 


NL 


NUM3ER 0" 
FLUORESCING 
ELEMENTS 


element number and absorption jump ratio of fluorescing 

elements 


1 

25 

0 


2 

20 

0 


3 

13 

1 

5 13.902 

4 

26 

0 


5 

14 

0 


6 

8 

•0 



1 

LAMBDA 

EC 

SIGMA 

OMEGA 

U 

25 

2,1020 

6,539 

3286.5 

,2600 

3.059 

20 

3,3590 

4,038 

2938,7 

.1201 

4,953 

13 

8,3400 

1,559 

2731,5 

,0195 

12.827 

26 

1,9370 

7.113 

3391,4 

.2925 

2,812 

14 

7,1260 

1,338 

2748, 7 

, 0273 

10,881 

8 

"0 

.327 

2694.7 

.0018 

24,198 


LAMBDA = WAVELENGTH OF EMISSION LINE 

Z * ATOMIC NUMBER 

R "s ABSORPTION JJMP RATIO 

EC * CRITICAL EXCITATION VOLTAGE (KV) 

OMEGA a K FLUORESCENCE YIEuD 

U 5 EO/EC 


Figure 12. Output from Program EMX2. Table show- 
ing what elements will fluoresce what other ele- 
ments and the calculated values of the parameters 
required for the corrections. 
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PROBLEM NO. 13 SPESSARTlNE RUTHERFORD PEGMATITE ARCH MAR 1966 20 K\ 


NIT 

2 

CON 

DEICDN 

r < chi > 

I/I(FU 

0 

25 

,278500 




0 

20 

, 007315- 




0 

13 

, 0 941 SIT 




0 

26 

,034820 




0 

14 

,176000 




0 

8 

,378214 




0 

TOT 

, 96902$ 




1 

25 

,277636 

-.080564 

.967478 

0 

1 

20 

, 007258 

-.000057 

,8815 03 

0 

1 

13 

,115107 

,020927 

.488548 

,004368 

1 

26 

,034043 

,000023 

, 973821 

0 

1 

14 

,175498 

- , 0 0 0 502 

.529520 

0 

1 

8 

,395990 




1 

TOT 

U 006330 




2 

25 

, 277947' 

. 000311 

.966396 

0 

2 

20 

, 007290 

.000032 

.877580 

0 

2 

13 

,116129 

,001023 

.484339 

,004176 

2 

26 

, 034844 

,000001 

.973792 

0 

2 

14 

,175743 

,000245 

.528780 

0 

2 

8 

,397283 




2 

TOT 

1,009235 




3 

25- 

,277972' 

,000025 

.966309 

0 

3 

20 

, 007292, 

,000002 

.877322 

0 

T 

“ IT 

,1IST9T 

.QUITO'S 7 

.4840^9 

,00*172 

3 

26 

,034844 

.000001 

.973775 

0 

3 

IV 

,175776 

.000033 

,528680 

0 

3 

8 

, 397389’ 




3 

TOT 

1,009470 





Figure 13. Output from Program EMX2, Correction 
of a single data point. Columns from left to 
right are the number of iterations, atomic number, 
concentration, the change in concentration from 
the previous estimate, the absorption correction 
factor and the fluorescence correction factor. 
Output from EMX gives the first estimate, and 
three successive iterations are shown. 
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PROBLEM no, 18 SPSSSARTINE RUTHERFORD PEGMATITE ARCH MAR 1966 20 KV 

CONICENfRATi ONS CSKHKTE'D FOR ABSORPTION AND FLUORESCENCE 


TOTAL 

25 

20 

13 

26 

14 

8 

93,3763 

18,7288 

1,1547 

10.9273 

8,7990 

16.6490 

37,1180 

96,1455 

18.5739 

1 . 1493 

11,1592 . 

10,1428 

17.0273 

38.0929 

96,5022 

18,6034 

1,0564 

11,3131 

10,3927 

16.9508 

38.1858 

97,8950 

18,1359 

1,2428 

11,4043 " 

11,2786 

17,1186 

38.6648 

97.1211 

13,5622 

1.3101 

11,2639 

10,6667 

16,9777 

38.3404 

99,9837 

16,0280 

1,5777 

11,4635 

14,0524 

17,4111 

39.3509 

98.8205 

18 . 4510 

1.2818 

11,3624 

11.4714 

17,2868 

38.9671 

99,6544 

19,3027 

1,1954 

11.5879 

11,1720 

17.1999 

39.1963 

99, 6049 

18,6356 

1.3306 

11,4349 

11,5640 

17,3604 

39.2295 

96.8116 

21.9537 

1.0463 

11.1276 

7,4876 

16.9910 

38.2055 

97,6590 

18,6423 

1.1641 

11,2189 

10 , 6568 

17,3221 

38.6548 

98,9539 

22,4654 

1,0229 

11,3598 

7,6522 

17,3949 

39 . 0588 

98.2569 

21,4451 

1,0600 

11,2451 

8,3259 

17,3567 

38.8241 

97.5696 

21 i 7993 

1,0403 

11,88114 

7,7861 

17,3046 

38.5588 

96,9564 

22,0603 

. 8808 

11,1623 

7,3568 

17,1506 

38.3456 

98,4979 

21,9154 

1,0781 

11,0927 

8,2432 

17,3580 

38 . 8105 

99.1520 

22,9757 

1.0419 

11.1640 

7,4123 

17,4843 

39.0736 

97,6577 

22.4470 

1,0264 

11,0772 

7,4208 

17.1865 

38 . 4998 

97.1199 

22,3219 

.9370 

10.9100 

7,3624 

17.2519 

38.3367 

98,2413 

25 . 4750 

.8093 

11 ,0863 

4,6991 

17.4026 

38.7691 

99,9353 

25,2997 

.8614 

11,6012 

5,1925 

17,5380 

39.4925 

100,8515 

26.5903 

.7972 

11,4786 

4,5605 

17.6944 

39.7305 

101.2307 

24,9433 

.9127 

11,6692 

6,0422 

17,7306 

39.9326 

TWi’ZWS 

2579*7 4 

.8111 

iii'sma ’ 

5,0323 

17.4663 

39.4503 

100.3386 

24,5256 

. 9095 

11.5202 

6,2801 

17,5563 

39 , 5468 

101.0406 

27,4685 

.7035 

11T6526 

3,6556 

17,7046 

39,8559 

100,9470 

27.7972 

. 7292 

11,6197 

3,4844 

17.5776 

39.7389 

98,9010 

26,6034 

.7529 

11,3837 

3,9621 

17,2473 

38,9516 

98,1012 

26,6759 

.7272 

ii’, 3295 

3,7342 

17.0314 

38.6030 

96,2212 

•26,9678 

. 6489 

11,0179 

2,8647 

16.8257 

37,8961 

96.6407 

26.6923 

.6313 

11,0732 

3,0250 

17.0544 

38.1645 

95.8910 

27.7789 

.5736 

10.9251 

2,0727 

16.7694 

37.7513 

95,4990 

25,9379 

. 6497 

11,0479 

3,1093 

16.9357 

37.8186 

94.7149 

27.2225 

,5648 

10,7778 

2,0112 

16,7402 

37.3903 

94,6782 

26.6446 

,5578 

10.6243 

2.8375 

16.7231 

37.2908 

95, 3665 

“ 27,4145 

.4743 

10,7855 

1,9065 

17,0524 

37.7333 

94.5974 

28,6018 

,4945 

10,7242 

,8421 

16.6573 

37.2775 

96.6447 

27,1306 

. 6296 

10,8672 

2,7010 

17,1696 

38.1466 

96,3126 

26T7688 

.6261 

10,8474 

3,2166 

16,9424 

37.9112 

96.1899 

28,4059 

. 4605 

10,8219 

1,2692 

17,1991 

38.0333 

97,7583 

26,4536 

. 6988 

11,3771 

3,6209 

17.0500 

38.5580 

97.4588 

28,5788 

.5460 

11.2490 

1,0957 

17,3591 

38,6303 

98.9293 

29,7418 

.4372 

11,4558 

,7668 

17,4308 

39,0969 

97.7847 

£9,0691 

. 5638 

11,3884 

,9550 

17,1653 

38.6431 

98.9956 

29,2985 

,5521 

11.5592 

,9335 

17,4626 

39.1896 


Figure 14. Output from Program EMX2. Concentrations in weight percent identi- 
fied by atomic number after correction for absorption and fluorescence. 
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PROBLEM NO, 18 S^ESSARTlNE RUTHERF ORD . PEGMAT I T E ARCH MAR 1966 20 KV 

concentrations coptrected tor ABSORPTION and fluorescence 


TOTAL 

MNO 

CAO 

AL203 

FEO 

S 1 02 

93.3768 

24,1908 

1.6157 

20,6482 

11,3199 

35.6122 

96,1455 

23,9308 

1 .6081 

21,0864 

13,0488 

36.4214 

96.5022 

24,0189 

1.4781 

21,3772 

13,3703 

36,2577 

97 . 8950 

23.4798 

1,7390 

21.5496 

14,5099 

36.6167 

97.1211 

23,9657 

1.8331 

21,2843 

13,7228 

36.3152 

99.8837 

20 ,6938 

2.2076 

21,6615 

18,0784 

37.2424 

98.9205 

23,8220 

1.7935 

21,4704 

14,7580 

36.9766 

99.6544 

24,9218 

1.6727 

21,8966 

14.3728 

36,7906 

99, 6049 

24,1249 

1,8617 

21,6073 

14,8770 

37.1338 

96,9116 

28,3444 

1.4640 

21.0267 

9,6328 

36.3437 

97,6590 

24 ,0690 

1.6288 

21,1992 

13,7100 

37.0520 

98.9539 

29,0051 

1.4312 

21,4654 

9,8445 

37,2076 

98,2569 

27,6978 

1.4832 

21.2487 

10,7112 

37.1261 

97,5696 

28. 1451 

1,4556 

20.9375 

10,0169 

37.0145 

96.95-64 

28.4321 

1.2325 

21 , 0922 

9,4645 

36,6851 

98.4979 

28,2950 

1.5084 

20,9608 

10,6049 

37,1288 

99.1520 

29 , 6640 

1.4578 

21,0*954 

9,5360 

37,3988 

97.6577 

29,9953 

1.4361 

20,9314 

9,5469 

36,7620 

97. 1199 

"23,8198 

T.31U 

*20 i &I55 

9,4717 

36.9017 

98.2413 

32,8907 

1,1324 

20.9486 

6,0453 

37,2242 

99.9953 

32,6645 

1.2053 

21,9216 

6,6802 

37,5138 

10H. 8515 

34.3307 

1.1155 

21 , 

5,8671 

37,8483 

101.2307 

32,2043 

1.2771 

22,0502 

7,7733 

37,9258 

~nm,22S5 

~ 33', 5235 

1,1350 

21 77128 

6,4740 

37.3603 

100,3386 

_ 31,6650 

1,2726 

21,7686 

8,0794 

37,5530 

101,0406 

35,4646 

.9843 

22,0187 

4,7029 

37.8702 

inn .9470 

35,838*9 

1 .(3263 

21 , 9565 

4,4827 

37,5986 

98,9010 

34,3477 

1,0534 

21.5107 

5,0972 

36,8920 

w;mr t 

~3~4 , 4 412 

“ rmirs 

' "21 , 4o63 

4,8041 

36,4301 

96,2212 

34,8181 

. 9080 

20,8195 

3,6854 

35,9902 

96,6407 

34,4625 

. 8833 

20.9239 

3,8916 

36,4794 

'95,9910 

'35.8654 

.8026 

20,6441 

2.6665 

35,9124 

95 . 4990 

33.4884 

, 9090 

20,8761 

4,0001 

36.2254 

94,7149 

35,1470 

.7903 

2n , 35*57 

2,5875 

35,8244 

94,6782 

34 ,4009 

.7605 

20,0757 

3,6505 

35,7707 

95 , 3665 

35,3949 

. 6636 

20.3802 

2,4527 

36.4751 

94,5974 

36,9277 

.6919 

20,2644 

1,0834 

35,6299 

96.6447 

35.0284 

.8810 

20.5347 

3,4749 

36.7258 

95.3126 

34,5512 

. 8761 

2n . 4973 

4,1381 

36.2399 

96.1899 

36,6749 

, 6443 

20,4490 

1,6328 

36,7889 

97,7583 

34.1542 

. 9778 

21,4982 

4,6583 

36.4699 

97.4539 

36,8981 

,7639 

21,2561 

1,4096 

37.1311 

98,9293 

38.3997 

.6117 

21,6470 

,9865 

37,2844 

97,7847 

37,5311 

.7889 

*21.5196 

1,2286 

36,7165 

98,9956 

37.8273 

. 7725 

21,8423 

1,2010 

37.3525 


Figure 15. Output from Program EMX2. Concentrations expressed as oxides 
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Figure 16. Output from Program EMX2. Plot of concentrations of oxides after correction. Similar plots 
can be obtained as the output of Program EMX. The analytical data in this problem are from Sinkankas and 
Reid (6) . 
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PROGRAM EMX 

COMMON SUMIT.NpROB, I D < 6 ) , NEL , MBKGO , NR , NO ( 1 3 ) * WTPC ( 10 ) , T AtJ ( 9 ) , 

1 NDM(0),TJMi^8OO) # CU^RBNt(8QrrrcouRrS(WI5,3),^R(adO>3), 

2 X INFO (80 0,8) ,KD(800,3>,NL<800,3>, I STOP , NUM { 26 ) , MNUM , NAME < 8 0 0 ) , 

3 NAmE2(50(5), kEY ( 9 ) > NP<9>, KEY2, IO'x, NS, IG, TOTAL(500), 

4 YMAX, 0(9), S I { 9 ) 

COMMON/i775/CONf5OQ,9),DlNFO(2OO.7),WD('20Q»3> 450 0 

EQUIVALENCE(TIM(?,XlNFO), ( COM, D INFO) , ( CON ( 1401) ,NLD> , ( NAME2, T IMS), 

1 ( TOTAL, CURRENT ) 

DEAD ( T , OC ) s OC/U.O - OC*T) 

NP AGE = 0 

JB = 22236420209020208 
C 

C READ AND AVERAGE D Af A FOR ONE PROBLEM 
C ' " 

1 CALL EMXINPUT (XMAX) 

I F { ISTOP) 95,2 ’ 

2 I F ( . NOT . NS ) 55 , 3 

3 CALL rmxhead 
PRINT 200 

200 FT5RT4A r C 25X ,13ffAvERAGED DATA/) 

•PRINT 201 

2 0 i ~F ORM A T ( 15X , 6H rEL . , 16X , 14HSPECT ROMPER 1, ‘ 

1 llX,14HSPECTRnMETER 2 , 11X , 14HSPECTR0METER 3/l5X,16H TIME CURR 
2ENT,3x,3(17HCtS/$EC SIGMA, 0X)7) 

C _ 

C FIND RELATIVE TIME, "PRINT AVERAGED DATA AND STAnUArD DEVIATIONS. 

C 

DO 5 M = 2,MNUM 
Ml = NUM(M-l) 

M2 ’= NUM(M) ““I 
T = 0 . 0 

DO 4 I = M 1 #M 2 
T = T + T I ME ( I ) 

Ti M5( I ) a T • . 5*T I ME ( I ) 

PRINT 205, I, NAME ( I ) , TIME{I),-CURRENT(I), (COUNTS! I, J), ERR ( I , J 
1 ), KD ( I , J ) , WL( I, J), J*i,3> 

205 FORMAT ( 1 5 , 2X , A8 , F6 . 0 , Flo , 4 , 3 < F10 . 1 , F10 . 1 » 13, 12)) 

ITIXMODFU, 50) 54*21 ‘ ‘ 

21 CALL EMXHEAD 
PRINT 200 
PRINT 201 

4 CONTINUE 

5 PRINT 310 
310 'FORMAT () 

I F < , NOT * NR ) 56 , 53 

53 I A= 0 

55 CALL EMXHEAD 
IFIRST = IA* I 
PRINT 203 
RRTNT 2D I 

203 FORMAT(l4HODEADf IME DATA) 

DO 54 I A= IFIRSf , NR 

PRINT 20 4, DINFq( IA),DINF0( I A , 2 ) , D T NFO (I A, 5 > , NLDU A ) , D I NFO ( I A, 3) , 

1 DINFO! IA,6),NLD( I A , 2 ) , D I NFO U A , 4 ) , D I NFO ( I A, 7) ,NLD( I A, 3) 

204 FORMAT(2lX,F10.4,3(2FlO f l, 15)) 

ITr;NoT.XMODFm>50} )55,51 

54 CONTINUE 
C 

C COMPUTE DEADTIME CORRECTION COEFFICIENTS WHERE NECESSARY AND CORRECT FOR 
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C DEADTIME 
C 

56 LAST = NUMTMNUM) - 1 
DO 7 i = 1,NEL 
ENCODE (5,3ao,Np(L> )NO(L) 

32 0 FOHMAK2X, 13) 

KEY(L) = 0 
IF<TAU(L) >7*6 

6 CALL DEAD? I ME ( L , XMAX) 

7 CONTINUE 

i f ( , not * ns ) i > i7 

17 DO 13 NN 2 1#NE|. 

JL AST = 3 

MM = NN 

IF(TAU(NN) .LT. 0.0)13*8 

8 IF (TAtJ(NN) .LT. 1,0)10*9 

9 MH = TAU(NN) 

10 TT 2 TAU(MM) 

DO 15 JX= IiJLAsT 
DO 15 I X= IiLAST 
IFTNLC IX# JX) ,NE ,NN) 15, 11 

11 JLAST 2 J,X 

COUNTS ( I X , JX ) 2 DEAD ( TT ,.C0UNTS( I X * JX ) ) 

15 CONTINUE 
13 CONTINUE 

CALL FMXHEAD 
PRINT 275 

275 FORMAT ( 25X, 34HDATA AFTER CORRECTION FOR DEADTIME/) 

PRINT 201 

PRINT 276., <13# £ COUNTS ( I 3 , J3 ) , ERR ( I 3 * J3 ) # J3 = 1 * 3 ) > I 3 = 1 , l AST ) 

276 FORMAT ( I5,26x#3<2FiQ .1,5X) ) 

C 

C FIND TWO SUCCESSIVE STANDARDS AND INTERPOLATE BETWEEN TO CORRECT FOR DRIFT 
C 

DO 35 MX 2 2,MNUm 
MIX 2 NUM(MX-l) 

M2X = NUM(MX) - 1 
DO 35 JJ 2 1.3 
DO 16 II = MIX * m2X 

IF(KD( II* JJ) .EQ.2, AND.NLCli, JJ) )18 i16 

16 CONTINUE 
GO TO 35 

18 12 » 11 M 

NL1 2 NL( U# JJ) ' 

20 DO 22 14= l2,M2y 

I r < KD ( 14 # JJ) .eg. 2. AND.NLC 14, JJ) . EQ . NL1 ) 24 , 22 
22 CONTINUE 
GO TO 35 
24 12 2 14 

CALL DRIFT (U, t2, JJ# IFLAG) 

K1 2 T2 * f“ 

IF(KD(K1* JJ) . NE, 2)28* 27 

27 11 = K 1 
GO TO 18 

28 IFUFLAG) 30i 32 
30 12 = K1 * 1 

DO ' TT) 20 
32 II = 12 
GO TO 18 
35 CONTINUE 



CALL EMXHEAD 
PRINT 277 

277 F0RMAT(25X,3iHDATA AFTER CORRECTION FOR DRIFT/! 

PRINT 201 

PRINT 276, (15* (COUNTS ( 1 5* J5 ) . ERR! l5*J5)#J5sl#3)*I5*l,LAST) 

C 

C SUBTRACT BACKGROUND FROM SAMPLES AND STANDARDS 
C 

CALL R*K~GD(lf 
CALL BKGD ( 3 ) 

37 CALL FMXHEAD 
PRINT 210 

210 Format ! 25x* 6 'qh'c n n c E' nTrat Ions c orrec t'ED"T0r' “d e a d f l me, d rT ft r and' b a 

ICKGROlJND// 21X,36HSAMPLE CORRECTED CONC 

1EN- rElATI VE/^X » 5iHE SAMPLE " CURRENT” CTS/5EC ‘TRATlON 
2ERROR/) 

DO 39 JP= 1,3' 

DO 38 IP= l,LASf 

38 ERR ( I P , JP ) = £RR( IP* J D ) /COUNtS( IP, JP) 

DO 39 IQ = 1*50Q 

39“Ctmno. JP5"*T'*0 ‘ 

DO 50 MY= •?* MnUm 
M' iY= NUM ! MY-1 ) 

M2 Y = NUM(MY) - 1 
DO 50 JY= 1*3 
DQ 40 NNN= i , NEi 

' DO 400 IY= Ml Y , Mg Y ' “ 

IF(KD( IY* JY> .EQ.2. AND.NLUY».Jy) . EQ . NNN ) 42, 40 0 
400 CONTINUE 
40 CONTINUE 
GO TO 50 

42 STD * 1 , 0/CQUnT S ( I Y * J Y ) 

SE = ERR^ rV,JY)i*2 
PRINT 310 

PRINT 213* I Y, NNN* JB* CURRENT ( I Y ) , COUNT S(Ty * JY ) * WTPC ( NNN ) , FRR ( I Y* JY) 
1 ) 

213 FORMAT (I5*I2,2x*A8,Fl0.4,F10.1*Fi0.5*FlD,3> 

DO 45 IZ = MIY^MSY 

IFCKDUZ* JY) . OR , NinZ ,’JY)7 n1 . NNN ) 45 , 43 

43 KN = KEY(NNN) - KEY(NNN) * 1 
IFCKN.LE. 500)44, 430 

430 PRINT 24 0 * NNN 

240 FORMAT < //55H DIMENSION EXCEEDED FURTHER DATA OMITTED FOR ElEMEN 

IT* 13//) 

GO TO 40 ‘ 

44 CON(KN*NNN) = CrUNTSUZ* JY ) *STD*WTPC < NNN > ' 

N AME2 ( KN ) = NAMFdZ) 

E = SQRTF(ERR.( 17* JY)**2 + SE) 

PRINT 213, l'Z, NnN , NAME ( IZ)*CURRFNT( I Z), COUNTS! 1 l * JY ) , CON ( KN * NNN ) , E 

45 CONTINUE 

TSO'TO '40 ~~ ' 

50 PRINT 310 

KEY! = KEY2 = KrY(I) 

DO 58 JZ= 2,NEL 

IF! KEY! JZ) . EG.Kf?Y2)58, 57 

57 PRINT 250, JZ 

250 FORMAT ( ////77ft NUMBER OF SAMPLE “DATA POINTS OF ELEMENT 1 DOES NOT 
1EGUAL NUMBER FOR ELEMENT *11) 

KEY2 = XMAX0F<Kf?Y2*KEY! JZ) ) 

58 CONTINUE 
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DO 62 I = 1, KEYS 
TOTAL ( I ) = 0.0 

DD”60 J = l,inr * “ 

60 TOTAL(I) = T 0 T A| 1 ( J ) * CON(I,J) 

62 'CONTINUE ’ " 

CALL PPLOT <NP> 

DO 72 I = i,KgY{? 

DO 71 J = 1,NEL 

ir<co\r(i;j) t LT.0'.o)6T,;i 
69 CON ( I , J) = 0.0 

71 CONTINUE ' 

72 WRITE (3*214) N *ME2 ( I ) , < CON ( I ,.J ) , Js 1 ,_NEL ) 

214 r OR M AT‘< A8 » 9 ( 5 pE~8 . 0 ) ) 

END FILE 3 

iruoxmri' 

91 DO 92 IC= 1,KEY2 
TOTAL ( IC ) = (f .15 
DO 92 JC= liNEL 

C O N ( I C > J C ) a CON( JC, JC)*(1,0 + O(JC) + Sl(JO) 

92 TOTAL ( I C) = TQTaLCIC) + CON(IC.JC) 

CALL PPLOT xnm) 

GO TO 1 
95 END FILE 3 
END 

SUBROUTINE EMXlNPUT (XMAX) 

COMMON SLIM!T,NPR0B, I D < 6 ) , NEL , HBKGD , NR , NO ( 1 0 > * WT PC ( 1 0 ) , TAIJt 9 ) . 

1 TyTDMr?) .TTRKbDTJ) , CURRENT < Bn U ) ,WUN T STB'0'0 T > EAR ( 8 0 0 > 3 ) , 

2 X I’NF 0(800*8) ,KD<800,3),NL<800,3), I S TOP , NUM ( 26 ) * MNUM * N AME ( 8 0 0 ) , 

3 NAME2C500)* kEY(9>* NP<9)* KEY2, tQX# NS, IGj TOTAL(BOG)# 

4 YMAX » 0.(9), ST (9) 

COMMDN/1775/D!NfO(200, 7) ,NLD(200* 3) , DUMMY ( 1995) *X( 101,5) 4500 

EQUIVALENCE ( T I ME , X I N^O ) , ( NAME2 , T I ME )#( TOT AL » CURRENT ) 

C 

C THIS SUBROUTINE REAqS DATA FOR ONg PROBLEM AND CALLS S.R. EMXSUM TO AVERAGE 
C 

ICHECK = 11UH11U112020B 

1 READ 100, NPRQB , ID, MEL * NDT., NS, MSKGDi SL I M I T # IQ* XMAX. YMAX 
100 FORMAT ( I*6»6A8»4.I2*F3,0, I2*2F3.0> 

I R RPR OB) 5, 3 " 

3 ISTOP = 1 
RETURN 

5 CALL EMXHEAD 
ISTOP = IOX = 0 

PRINT 200 , REL » NDT, NS, MB KGp, SU MlT 
2DirrOWWT' riBFTW.' DF"'ECEtfENT5 8 , 1 2/ 1 3,.26H DEAD! TR'E* “RES D I7QT3S PFFT SET / 

1 1 3 , 24N SAMPLE READINGS PER SET/28H BACKGROUND READING METHOD =, 

2I2/39H SIGMA FOR STANDARD PEAK MAY* NOT 'EXCEED# F5.1#11H * SGRTF(N) 

3 ///10X, 50HATOMIC PER CENT W T AU OXYGEN SILICON/ 

4 60 H ELEMENT NUMBER "OF STANDARD FACTOR ‘FACTOR 

5 /) 

read - - rtr2, morn * wtpctt t tttottttot rr mrn » nomt p » i * nel ) 

102 FORMAT (IX, 1 3* F6 , 0 4 E15 . 5-, 2F10 . 0 * 5X . A5 ) 

PRINT 20 2, ( I,NOU ) * WTPC( I >TTSUm * OT I ) , SI XT) » NOM ( I ), 1 = 1# NED 
202 FORMAT < 15, 1 9 , 5x * 2PF10 . 6# OPEli . 2 , 2FlO » 5 » 3X * AS ) 

DO 907 I = 1#NEl 
IF ( 0 ( I ) >908, 9Q7 

9B7 "CUNTtNUE 

GO TO 1111 
908 IOX' * 1 

DO 909 I = 1#NE{. 



lF(SI<I>>inOO,9fl9 
909 CONTINUE 
GO TO Uil 
1000 IOX b - 1 OX 
1111 CONTINUE 

CALL EMXHPAD 
PRINT 204 

204 F0RHAT(23X»UH JNPUT DAJA/5X,68H _NAME TIME 

1SPEC .1 SPEC 2 SPEC 3” KOBE/ ) 

N = NR = 0 

IB = 12121212121220208 " 

M = NUM(l) « 1 

7 N = N + i 

IF(N.GT , 80 0 ) 9 0 ^ A 

8 READ ( 2,'l04)‘"(X<l # j) , Jsl, 5 jYkODE, N AME ( N ) 

104 FORMAT (F6.2/F6.3 i 3F6.0, A6,A8) 

9 !F(KODE.EQ. iCHEeK)96#10 

10 IF(KODE.NE. !B)ijni 

11 M = M + 1 
PRINT 310 

310" FORMAT ( )""" 

NUM(M) = N 
IF ( M * E Q • 2 5 ) 9 6 i 8 

12 PRINT 206, N, NaME(N), ( X ( 1 , J ) , J = 1 , 5 ) , KODE 
206 FORMAT ( 1 5 » 2 X , Aft , pl 0 . 2 » F i 0 * 3 » 3F 1 6 , *0 ^ ?X , Aft ) 

DECODE (B,3 Q0,Kf) D5) <K D(N,J) , N L ( N , J ) * Jal , 3 ) 

300' FORMAT (611, 2X) 

IF ( KD( N) .GE, 2)25, 13 

13 IF( NS, NE. 1)17,14 

14 XINFO(N) = X 
X ■= 1,0/x 

DO 15 J = n 3 

K = J + 2 

15 ERR ( N, J ) = SQftTF ( X ( 1 , K) )_»X 

DO 16 J = 2,5 

16 X I NFO ( N , J ) * Xd,J>*X 
GO TO 7 

17 , D 2_?° 1 3 ?>NS 

READ (2,i04)" ( X (~I i J )*, J= 1 , 5 ) rJKODE, JUNK 
PRINT 206, N, JUNK, < X ( ! , J ) , J = 1 , 5 ) , JKODE 
IF( JKODE, EQ,KQDp>20, 19 

19 PRINT 208 

208 FORMAT (/ 75 1H ERROR \U GROUPING DaTA--Fu'RTHE'R "PROCESS 
65 IF( JKODE. EQ. ICHPCK)1, 95 

20 CONTINUE ~ “ 

CALL FMXSUM (N> NS> 

GO TO 7 

25 IF(KD(N) .NE, 4)20,-26 

26 NR = NR * 1 
IF(NR,GT. 200)90, 27 

27 do 50“ r s t>;'ndt ~~ 

READ (2,104) ( X ( I , J ) , J-l, 5 ) , JKODE* JUNK 
PRINT 206, NR,JuNK, C X (I , J ) , Jsl , 5 V, JKO'DE 
IF( JKODE. NE.K0DPU9, 50 
50 CONTINUE 

CALL EMXSUM (N*NDT> 

DO 55 J = 1 ,7“ ’ 

K = J + 1 

55 D I NFO ( NR , J ) = X f NFO ( N , K ) 

DO 60 J = 1,3 


CURRENT, 
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60 NLD(NR,J) = NL(NtJ) 

GO TO 8 

28 DO 61 J = 1,3 

I F< KD( N, J) . EQ . 2 ) 62 , 61 

61 CONTINUE 
GO TO 13 

62 I = 1 

29 I = I + t 

I F ( I ,GT,iOi)90^3O 

30 READ (2,104) ( X ( I , J ) , J='l, 5 ) , JKODE, JUNK 
IF( JKODE, NE,KODj?)31»75 

75 PRINT 206, N» JyNK , ( X ( I , J ) , J = 1 , 5 ) , JKODE 
GO TO 29 
35 NX * I - 1 

CALL EMXSUM ( N* NX) 

1 F ( ,NDT. IST0P)3p, 64 
64 ISTOP =0 
GO TO 65 
32 KODE = JKODE 
N = N + 1 
IFTN.GT. aDD)9T5L42 

42 NAME ( N ) = JUNK 
DO 43 J = 1,5 

43 X< 1. J) = X ( I , J) 

GO TO 9 

90 PRINT 250 

25 G FORMAT (///52H ERROR- -DT MENS I ON FXCFEBED- -DATA CANNOT BE PROCESSED 
1/) 

95 READ 150, JKODE 
150 FORMAl(30X, A6) 

IF(JK0DE,E0, I CHpCK ) 1 * 95 

96 MNUM = M + 1 
NUWMNUM) = N 
END 

SUBROUTINE EMXHgAD 

COMMON SUMIT,NPR0B, I D ( 6 ) , NEL , MBKGD , NR , NO ( 1 0 ) # WTPCC 10) , DCA( 9) , 

1 NOM( 9) ,TIME(flOO)# CURRENT (BOO), COUNTS ( BOO »-3), ERR (800*3), 

2 XINFO(800,8),KD(800,3),NL(800,3), I STOP , NUM ( 26 ) *# MNUM , NAME ( 80 0 ) , 

3 NAME 2 150 0 ) , REYC9), NP(9), KEY?, IOX, NS, IG, TOTAL f 50 0 ) , 

4 YMAX, 0(9), S M 9) 

EQUIVALENCE ( T I ME , X I N r Q ) , ( NAME2, TIME), ( TOT AL , CURRENT ) 

C 

EV> C PAGE HEADING ROUTINE 

ro c 

4^ 'N'P'AGE = NPAGE * 1 

LK PRINT 200, NPROR, ID* NP AGE 

H- 200 F0RMAT(12H1PR0B|.EM NO. I 6 , 5X , 5 A8 , 30 X , 4HPAGE 13/) 

END 

SUBROUTINE EMXSijM (NY, NX) 

COMMON SL I M I T , NpRQB , I D ( 6 ) , NEL , MBKGD , NR , NO ( 1 0 > * WTP C ( 1 0 ) , DC At 9) . 

1 N0M'C9),TIM^t BOO), CURRENT! 800T;cnUNTS(«00, 3), ERRT 80 0,3), 

2 XINFO<800,8),KD(800,3>,NI(800,3),I STOP, NUM ( 26) , M NUM , NAME ( 8 0 0 ) > 

3 N AME2 (500 ) , kEY<9), NP(9), KEY2, IOX, NS, IG, TOT AL t 50 0 ) , 

4 YMAX, 0(9), $1(9) 

COMMON/1 7 75/DUMMY ( 399 0 ) , SUM (5), X( 101,5) 4500 

EQUIVALENCE ( J I ME , X I H? 0 ) , ( NAME2/T I ME ) , ( TO! AL ^CURRENT ) 

THIS SUBROUTINE AVERAGES DATA SETS, COMPUTES S 1 ANDARD DEVIATIONS, 

AND CONVERTS TO CQUnTS PER SECOND. 
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N = NY 

XN = FLO ATF ( NX ) 

.DO 5 JJ= 1,5 
SUM C J J ) = 0.0 
DO 5 I = 1, NX 

5 SUM(JJ) = SUM ( J j ) + XU,JJ) 

T I ME ( N ) = SUM 

X = 1. 0/SUM 

CURRENT! N) s S U M ( 2 ) * X 

DO 20 J = 1,3 

K = J + 2 

IF ( NL ( N, J ) ) 9, 7 

7 counts ( Nr j) = err('n, j) = d.o 

GO TO 20 

9 TEST, s S Q R I F ( 5 U M ( K > ) * X " 

S = SUM ( K ) / XN 
DEV = 0.0 
DO 10 I = 1 , NX 

10 DEV = DEV + ( X ( T , K ) - S)**2 
DEV = S0RTF(D6VJ*X 

■ " TFTKD ( N ; JJ7NE721 15,11 

11 IF(DEV.LE.SLIMIT*TBST)15i 12 
1? DEV = DEV/TE'Sf 

PRINT 200, DEV 

200 FORMAT! ////32H SCATTER fN STANDARD READING IS ,F4.1,33H * SGRIF(N) 

1- - FURTHER DATA OMITTED////*) 

I STOP' = '1 
RETURN 

15 ERR*(N,J) = MAXlFt DEV, TEST) 

X INFO { N , K ) = SUm(K)*X 
20 CONT I NUE 

END _ 

SUBROUTINE DEAD? I ME (N, XMAX) 

COMMON SLlMl*T,NpROB, I D ( 6) , NEL, MBKGD# NR, NO ( 10 ) * WTPC ( 10*> , T AU < 9 > . 

1 N0MC9) , TIME! 800 )•, CURRENT ( 8n0 >, COUNTS! 800, 3), ERR ( 800, 3) , 

2 X I NF 0T8 0 0 » 8 ) » K D ( 8 0 0 , 3 ) » NL(80O>-3) , ISTOP , NU» ( 26 > , MNUM , NAME ! 80 0 ) , 

3 NAME2!5(J0), kEY!9), NP(‘9->. KEY2* IOX, NS, IG, TOTAL<500)» 

YMAX, 0 ! 9 ) , SK9 ) 

’ C0MM0N/177 5/D I NFO ( 20 0, 7) , NL 0( 20(1 , 3 ) ,.DELY1 20 0 ) , Y( 20 0 ) , B ( 2 ) , SB! 2 ) , 

1 DUMMY (2096) 4500 

EQUIVALENCE (dCfS,DlNt'O(20i) ) , (b'E,fiiNFfl<801 >> 

EQUIVALENCE ( T! ME , X I N- 0 ) , ( NAMS2 , T I ME ) , ( TO TAL , CURRENT ) 

DIMENSION DC Ts ( ?0 0 , 3 ) » D"E (.20 0 , 3 ) 

DO 5 J = 1,3 

W "5 U =‘i,NR “ 

I F ( NLD ( II, ‘J) , EQ . N ) 6, 5 
5 CONTINUE 
PRINT 250, N 

250 FORMAT !.///?9H No DEADTIM"E DATA FOR ELEMENT , \ 2 » 38H NO DEADTIME CO 

1RRECT I ON WILL Bp MADE) 

T'AXTTN) = -r,0 
RETURN 

612= T 1 + 1 

DO 8 I = 12, Nr 
IF(NLD(I,J-).EQ.N)8,10 

8 CONTINUE 

12 S NR 

ND = NR ~ II + 1 
GO TO 11 
10 12 = I - 1 
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ND = I - 11 
11 DO 15 I a U, 12 

~ -UCTS1 I , s ( I , JT7DTNTD7T7 " ' 

15 DE ( I » J )* = ( D 1 NFo < I )-/DH ( I # J ) ) **2 

CALL EMXHE A'D ‘ ' ' 

CALL l.SO (ND#i#SlG^DlMFOUl)#DCTS(Il#J),DE(Ii#J)»YC.Jl). DSLY,R,S8> 

TAU(N) = -B(2 >71*** 

PRINT 200, NO ( N ) » J, TAU(N) 

2T0 FOTR*f \ 777Z7T5Wl?m ’EL'EmENT" HETTER ON SPECtRITmITIr, 12* 

1 8H, TAU = ,E9«.2r8H SECONDS ) 

DX = XMAX* , 10 ' *■ 

DY = ,125*8 

D I V = IDOL * XMaX ” ‘ 

CALL CCAXIS ( 0-* , , 14HC0UNTS/CURRENT# 14> 8 *,9O*f0.»DY>12,5) 

uao: ■ ruA'xrsc o ;;"20 -rouititent' r^rcfto - amps j , - 20 » rirv ,r. , 5/. ;w;-tnvr 

call CCAXIS, , -1 , 8 . , 90 . „0 . , Dy , *12 f 5 ) 

CALL CCAXIS 0 t ,20HClTRWEWT ( M lWO-WPST'* -20V10 . , 0 . , 0 . , DX > D I V J 
CALL CCS YM80L ( 1 . * 10 , , . 21 , I D , 0 . , 48 ) 

CALL P2LINE CD InFO ( 1 1 ) , D'CTS ( M , J ) , NO , 1 , ->1 , 3, . 07 * 0 . 0-, DX , 0 . 0 , DY ) 

YNEW = (B + XMAy*9(2) )/DY 

C“ALX ' c C POFTTUr; , YNEW 7^7 

CALL CCPLOT < 0- , 8., , 2 > 

CALL CCPLOT 

END 

SUBROUTINE LSQ ( M , I W * S I GM A , X , F2 , W , Y , tTELY ,.B , SB ) 

C 

C STTOHT VERS-IOM OF F2 UCSD LWOL — FOR FI TTTNE HTTH STRAIGHT LINE ONLY, 
C F TN- 63 
C JANUARY. 1966 
C 

DIMENSION 5(2), X(2Q0), F2 ( 20 Q ) , *ST<2), SBC?)*. F C 20 0 > > ’P M(200 W 
1 ‘P < 2 0 0 ) , 8(2), DELYC 200 ), W( 200 >, A<2,2)j T(2), Y( 200 ) 

LL = 0 ' * “ 

A = A (2, 2) = 

FBAR 3 xBAP = F m = 0,0 

DO 10 I = l',M 

IF(IW)1009,i010 

1010 W2 = W(I ) = 1,0 

GO ' TO* 1011 ' " ’ " 

1009 W2 = SQRTFCW ( I ) ) 

1011 FM = FM + W( I ) 

F { I > = W2*F2< I ) 

PO PM ( I > = W2 

! Nj FBAR 3 FBAR + F< I >*PM( I ) 

■£* io KB**-"* -xbar -t xrn'»ptirrTT*2 ■ — ■ — ~ 

LN X8AR 3 XB’AR/FM 

LH T(l) 3 FBAR/FM 

A( 2, 1 ) = -X8AR 
PXF « PXP 3 0,0 
DO 20 I 3 1 ,.M 

-PCI) = t X < n—yBARI^PMTiT —■ ■ 

PXF 3 PXF + P( I S*F( ! ) 

20 PXP » PXP + PCI )**^ 

T ( 2 ) = PXF/PXP 

SCI) = PMXPM = FM 

BCD = TCI) * A ( £ » 1 ) * T ( 2 ) ★ A ( 2 , 1 ) 

- Tt?-) m2,“2i “ 

165 SIG2 3 0,0 

DO 180 I = 1» M 

YC I ) 3 8(1) * B|2)*X< I ) 
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175 DEL Y C I > * Y(I> * F2( I > 

180 SIG2 = SIG2 + (DELY(t)**2)*W(I> 

SIG2 = SIG2/(H - 2) 

SIGMA s SQRTF ( S t G2 ) 

S ( 2 ) = RXP 
DO 499 I = 1,2 

499 STM) = SIGMA/SrRTF(S( I) ) 

DO 501 I = 1,2 

SB ( I > = 0.0 “ 

DO 500 J =; 1,2 

500 SB ( I ) = SB(I) + ( A( J# I )*STC J) )**2 

501 SB ( I ) = SQRTF(SB( I > ) 

192 PRINT 1,(1, 1=1,2) 

1 FORMAT (36H0CoEFFiCI,E\|TS OF YaBl+B2X AND ERRORS// 2(3H B(I2# 

1 2H ) =El5 . 7 , 6H ERR8 S E10 • 3") ) 

185 PRINT 186, SIGMA 

186 FORMAT ( 8 H 0 S I G M A *E16.7) 

670 PRINT 2, ( 1 , X{ I ) , F2 ( I ) , Y( I ) , DELY< I) , W ( I ) # I S 1 * M ) 

2 FORMAT ( 47H0 i XM> FI ! ) 

137HY ( I ) DELY( I ) W( I )// ( I 6 , 5E16 . 7 ) ) 

190 CONTINUE 

END 

SUBROUTINE BRfFf U 1 #1 2 , J , I FLA G ) 

COMMON SLIMIT,NPR0B, ID<6> , NFL , MBKGD , NR , NO ( 1 0 ) # WTPC ( 1 0 ) , DC A ( 9 > * 

1 N0M< 9), TIM6’( 800), CURRENT (800) , COUNT'S ('80 0 ,3‘)V£RR ( 80 0 * '3 ) . 

2 XINFO(BOO#8> .KD(800,3) ,NL(a00,3> ,1STOP#N‘UM(26) » MNUM > NAME ( 800 ), 

3 NAME2< 500 ) , KEY ( 9? i NP{9), KEY2, IQX, NS," IG, f OT AL ( 50 0 ) > 

4 YMAX, 0(9), fi I ( 9 ) 

EQUIVALENCE ( T I ME , X I N r O ) , ( NAME2 , T I MS >,( TO T AL * CURRENT ) 

IF (ABSF( COUNTS ( fi. J) - COUNTS ( I 2, J ) ) , LE . MAX1F ( ERR ( U, J ) , ERR M 2, J ) 
1 ' >>90*2 

2 T = TIMEU1) 

CF = (COUNTS M2*. J ) /COUNTS (1 1‘, j)~- 1 . 0 )7 ( T I ME <T 2 ) - T) 

5 K1 s II + 1 
K 2 = J2-1 
I FLAG = 0 
DO '10 1 = Ki,K2 

IP COUNTS ( I * J ) = COUNTS! I ,J)/(1 . 0 + CF* < TJME< I >-T ) > 

RETURN ‘ ’ ‘ 

90 PRINT 200, I 1, I?, J 

200 FORMAT ( //31N DRIFT BETWEEN STANDARDS LINES ,13, 5H AND ,13. 

i 14 h, Spectrometer# 12 , 

1 51H IS WITHIN COUNTING FRROR--NO DRIFT CORRFCT*! ON MADE) 

IFLAG s 1 

END" 

SUBROUTINE BKQD ( KKD ) 

COMMON SLIMIT,NpROB, I 0? 6> ,NEL, MBkGU, NR~, NU( 10 )'» WTPC (10 ) , DC A (91 . 

1 NOM(9) ,TIME(800 ) , CURRENT ( 8 0 0 ) ,COUNTS(800,3),ERR(800,3), 

2 XINFO(8GO#8) ,KD( 80 0 , 3 5 , NL ( 80 0 , 3 ) , ISTOR , 26 ) # MNUM # NAME ( 80 0 ) , 

3 NAME2 ( 500 ), KEY ( 9 ) # NP(9), KEY2, IOX, NS, IG, T0TAL(5OO># 

4 ymax rum , s n 97 ** - — 

EQUIVALENCE ( Tl ME, X INFO) , ( NAME2, T I ME ) , ( TO T AL * CURRENT ) 

C 

C THIS SUBROUTINE LOCATES BACKGROUND READINGS AND SUBTRACTS THEM FROM 
C THE APPROPRIATE PEAK. 

C 

KOTTE- * KKD 

MODE = KODE - 1 
DO 50 J = 1,3 
DO 50 M = 2 , MNuM 
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Ml = NUM(M-l) 

M2 = NUM(M) - 1 

DTTmI ~ ~ 

3 DO 5 II = Ml, M2 

r f < kd ( ii, j) .eq.'kode. amd.nl ( u, J>Ye,5 

5 CONTINUE 
GO TO 50 

8 12 = UJL-I __ 

J2 ~=~ M2 

Nil a NU(IdjJ) 

DO 10 I = 12, M2 

1 F ( KO( I f J) .EQ.KRD6. AND.NK I , J) . EQ. NLl ) 12, 1 0 
10 CONTINUE' ' ' '■ 

b = counts ( ii, j > 

E B SG ' s "Ef? kill 7JY* * 2 
GO TO 17 

12 12 * I 

J2 = 12 + 1 
DO 13 I a J2, M2 

IF(KD( I , J) .EQ.KQUE. AND.NK I , J) . EQ . NLl > 15 . 13 

13 ~CONT77\IUF“" " ' ” 

J2 = M2 

GO TO 22 
15 J2 a I - 1 

22 B = ( COUNTS C tl-.J) + COUNTSU 2 , J ) ) * .'5 

EBSQ = ( ERR ( 1 1 , J ) **2 + ERR ( I 2 , J )**?)*, 25 

17 IF( JTOT.TOKfflJ73f.23 " ” ' 

23 DO 25 I = J1,J2 

IF(KD( I, J) .NE,MqDB.OR.NL< I, J) . NE . NL l ) 25 , 24 

24 COUNTS ( I , J) a COUNTSU, J) - B 
ERRU,J) = SORTf(ERRU, J)**2 + EBSO) 

25 CONTINUE 

Ji = J2 + 2 
GO TO 41 

30 DO 35 I a J1,J2 

IFCKDt I f J) - EG. MODE. AND. NL( 1 4 J) • EQ • NLl ) 40 # 35 
35 CONTINUE 

PRINT 205, II, f 2 

205 R3RWAT(/7/3^H‘TIo SAMPLE OAT'A FOR“B~A C KG ROUND', LINES ,13, 16) 

GO TO 41 

40 COUNTSU, J) = CflUNTSU, J) - B 
ERR(i»J) = SQRTf ( ERR ( I , J ) **2 + E8S0) 

J1 s I + 1 

41 I F < J2.GE.M2)50,42 

42‘ * 1' ■* ‘ 

GO TO 8 
50 CONTINUE 
END 

SUBROUTINE PP LOT (N) 

COMMON SLIMIT,NpROB, ID(6),NEL,M8KGO,NR,NO(10)»WTPCC10-),TAU(9) , 

1 NGMX 9) , TTMHT 80O’> i'CURRENTT8U 0 r,TTyUNT5T'8IJ 0,3), ERR ( 800,3), 

2 XINFO<800,8) ,KD<800,3) ,NL(80n,3) , IST0P,NUM(26) , MNUM , NAME < 80 0 ) , 

3 NAME2C500), KEY<9>, NP<9), KEY2, IOX, NS, IG, TOTAK5Q0), 

4 YMAX, 0(9), SI<9) 

COMMDN/1775/CON«500,9) , D I NFO <200 , 7) , NLDCZO 0,3) 

DIMENSION JPREF<9), X(500), N(20) 

EffUTV'A L END E T TT Ml ; XI NFT5 5',' TCDNTUTnTo ) , ( CON (147)13 , NlD 5 , (NAME2 , T f ME) , 
1 (TOTAL, CURRENT), (X, ERROR) 

DATA (IPREF = 4 ;3 , 5 , 0 , 1 , 2 , 6, 7 , 9 ) 

CALL FMXHEAD 
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PRINT 220 , (N(Ll,m*NEL> 

220 F0RMAT(25X»6QHCftNCEN-T RATIONS CORRECTED FOR DEADTIME* DRIFT* 
iCKGR0u‘ND//l5X, i 5HT0TAL»lf)(5x; AS ) ) ' ' 

DO 10 I = 1 * KEY? 

PRINT 256, NAME?( j ) * TOT A I, ( I > * ( CON< I > J J ) * JJ = 1*NEL ) 

256 F0RMATC2X,A8.11<2PF1Q.4)> 

I F C XMODF ( I .5) )1Q,5 
5 I F ( XMODF < 1,45) >3,7 

7 CALL EMXH6AD 

PRINT 220, ( N( K ) * K S 1 » NEL) 

GO TO 10 

8 PRINT 400 
400 FORMAT {) 

10 CONTINUE 

14 I F ( I G ) 1 5 * 9 0 
90 RETURN 

15 I F( N . NE . NP ) 30 * 16 ’ " ' 

16 IF ( l OX ) 90 , ,30 
30 YLEN = 8*6614 

YYMAX = YMAX*1 » i 

FLaXTFTTj<FY‘2+7y78 _ I ~~ 

PX1 = XLEN * l.Q 

PX2 * = XLEN + i . <5 ‘ - 

CALL CCPUO r (X lEn + 3, 0* 0, o > 

33 DO 35 I A - 1,KEY| 

35 X< I A) = I A 

PTl *."8924 “ ‘ 

PY2 = .7874 
DX = 8.0 
DY = YYMAX/YLEN 
DO TO~*J s l.TJEC"”’ 

CALL P2Ll'Nb_CX* CONU.J). KE Y ( J ) , 1 , I P_REF ( J-) , . 0 7 , 0 . 0 , DX . 0 
CALL CCS YMB'OL ( pXi * P Y 1 , . 21 , I PREF ( J ) , 0 . 0 , -1 ) 

CALL CCSYMBQL_ CpX2,PY2, , 2l*N(J)*0.0,5) 

PY1 = PYi + .5 

PY2 * PY2 + .5 

40 CONTINUE 

CALL P2LINF (X* TOTAL, KEY2, 1 , 1 , 11 , . 0 7 , 0 , 0 , DX ,0 , 0 , D Y ) 

"CAlL CCSYMHOL* (pXr.pVi * .'35*1170.0,-1) 

CALL CCSYMBQL (pX2,PY2, . 35 , 5HT OT AL » 0 . 0 * 5 ) 

*call ccaxis ( o • , o . » 60 hconcenTrati ons "Corrected for deadtime 
l, and background, 60, ylen* 90 . o, o . n , dy* 12. 7) 

CALL CCAXIS (XLpN, 0 . * ID, 8, YLEN, 90 . * O', *DY* -12. 7) 

CALL CCAXIS ( 0. ,7,874, ID*-48*XLEN, 0 .0* 0.0*8. 0 ,-8.0) 

CALL C*C AXT S FO'.Ttn * I'D', -48," XLEN, 0 . 0, 0-, 0, 8, 0 *8.0) 

CALL CCAXIS (0.,YLEN , ID* -48* XLEN, 0 . 0* 8 . 0# 8.0 ,-8.0) 

DO 50 IC= 1.KEY2 
IF(NAME2< IC> .E0.1H )50,45 
45 PX = Xi IC)*,125 

CALL CCSYMBOL (pX* YLEN* . 14* NAME 2 ( I C ) * 45 . #8 ) 

50 -CtmTTNUF ‘ ~ 

CALL CCPLOT ( 0, 0,-3) 

END 

SUBROUTINE P2lInE(X,Y,N,K,J,L*H,XMIN*DX,YMIN*DY) 

CLINE VERSION 5 CdPyRIGHT 1964 CALIFORNIA COMPUTER PRODUCTS 
DIMENSION X(1 ) , Y { 1 ) 

X is The NAME of the array of ordinate values. (SEE 
Y is The name of the array of abscissa values. (SEE 
N is the number of poi-nts in the array 


AND BA 


0 . DY) 


DRIFT 


LlNEOOiO 

L I NE 0 0 3 0 
NOTE) LINE0040 
NOTE) LINE0050 
L I NE 0 0 60 
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K IS TwE REPEAT CYCLE OF “A ‘ “HIXElT ARRAY (NORMALLY = i). LINE0070 
J IS THE ALTERNATE NUMBER OF DATA POINT TO PLOT A SYMBOL . L I NEO 080 
‘ J-WfLI = (TTOIFTTNE PLOT* NEQAV I'VE FU* 'POINT VC5T , ' lTN'EOO'Sq 
J = 1 FOR POINT FOR EVERY DATA POINTS FOR EVERY OTHERL I NE 0 10 0 
L IS Anj INTEGER DESCRIBING SYMBOL" TO Bg USED, SEE SYMBOL LlNEOliO 
ROUTINE FOR LIST LINE012Q 

LINE01&0 

NP = N * K * 1 LINE0170 

mr*lP‘*T ' ““ ‘ LINE 0180 

NO = NP - K LINE0190 

13 = 3 ‘ ‘ L I NE024 f) 

CALL CCWMERE(XN, YN) LINE0250 

12 = -1 * ’ ‘ LINE 0280 

XN = XN * DX * xM I N LINE0270 

YN 5 YN' * DY * 9MTN ' ' ‘ T. TNED2B0 

DX1 = ABSF ( X ( 1 ) *XN ) LINE0290 

DY1 = A8SFI Y(i )*TN> LINE0300 

DX2 = ABSF < X ( NO ) - XN ) LINEQ31G 

DY2 « ABSF ( Y ( NO ) - YN ) LIN6032O 

DX 1 = MAXlF'(DXl.DYl) UNE0330 

' U‘X'2~ 5 ~MA'X 1 F ( DX'2‘;TTY2 ) * LTNE0340 

KK = K l INE0350 

NA=-NT = J $ IF(N*f) 100,101,102 


100 NT=-NT 

102 NA = ? N-l ) /NT 

101 NA = NA * NT * NT +1 - N 

iriDXl-DZ?^ ID, 10, 20 ” - ---- 

10 NO * 1 
KK = -KK 
NA = NT 
20 NV = 2 
NW = -2 
KL = 2 

IF (J) 30,40,50 

30 NV = -3 

NW = -1 
KL = i 
GO TO 50 

AD W='NP * ‘ 

50 DO 90 I = 1#N 

IF (NA-NT) 65,80,70 

60 CALL CCSYMBOLI (X< NO)-XMI N) /DX, ( Y( NO) "YMIN) /DY,H#L, 0 . 0, 12) 
NA = 1 
GO TO 80 

85' 'WTO asmTTKl 

70 CALL CCPLOTMX(NO)-XMIN) /DX, ( Y ( NO ) - YM I N ) /D Y , i 3) 


75 

na = 

NA * 

1 

80 

13 = 

NV 



12 = 

NW 


90 

NO = 

NO - 

KK 


-RETURN' 

END 


LINE0370 

L INE03S0 
LINE039O 
LINE0400 
L INE0410 
LINE0420 
ll NE0430 
UNE0440 
L I NE 0 450 
l INE0460 
LINE0470 
L I NEO 480 
l INF0490 
L INE0500 
L I NE 051 0 
L INE0520 
LINEQ530 
LINE0540 
L I NEO 55 0 
L 1 NE 0560 
L INE0570 
LINE0580 
LINE0590 
L I NEO 6 0 0 
L INE0610 
L I NE 0 620 
LHVE0530 
L I NE 0 640 
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PROGRAM EMX2 

COMMON NPROB* I D ( 6 ) ♦ NEL*_EO* NZ<20), SCON(20), SCF(20), MUBAR* 

1 0(20)* IOX, NELP1* NF ( 20 ) * NZF<20*4), Z(20)* R(20*4), EC(20)* 

2 SIGMA( 20 ) 9 WK(20J* MU L2.(U2iU * JJl (-2 0 * 4 )j 51(20)-* NVAR* NELP2 » 

3 OCON ( 20 ) * C0N(20)> DELCON(20)* NOM(20)* NS* NP(20)* IMIN* 

4 IMAX* IG* YMAX, DMXN 

COMMON /PP/NA V E 2 ( 500 ) * Y(500*20)* TOT (500) 

NTOTAL = 5H TOT 
CALL CONSTANT 

1 NPAGE = IOX = 0 - 

CALL SETUP 

I F ( SCF ) 9 * 4 

4 DO 5 I = 1 * NEL 

5 CON ( I ) = SCON ( I ) 

CALL OXYGEN 
PRINT 200 

200 FORMAT ( 1 IX * 53H STANDARD FOR ELEMENT Z OXYGEN SILICO 

1N/94H NL Z CONCENTRATION F (CHI ) * ( 1 + FL IJOR ) FACTOR FAC 

2T0R F(CHI) I ( FL ) / I / ) 

DO 6 I = 1 > NEL 
FCH I = ABSORB ( I ) 

FL = FLUOR(I) 

SCF ( I ) = FCHI*( 1.0+FL) 

PUNCH 302 * I* NZ ( I ) * SCON! I)* SCF(I)* 0(1)* SI(I)*NOM(I) 

302 FORMAT (2I5»4Fl0.6,5X*A5 ) 

6 PRINT 202* I* NZ ( I ) * SCON ( I ) * SCF(I)* 0(1)- -HI), FCHI* FL 
202 FORMAT! 13* I 5 * F 13. 6 *F2u.6 * F 13 .6 * FlO . 6 * IOX* 2F10.6) 

GO TO 1 
9 NS = 0 

8 READ (3.100) NAME . (OCON ( I ) ♦ I = ]*NEL> 

]00 FORMAT (A8,9E8.0/8X,9E8.ri) 

I F ( EOF *3)75*7 

7 DO 11 I = 1 » NVAR 
CON ( I ) = OCON ( I ) 

11 OCON ( I ) = OCON ( I ) *SCF ( I ) 

CALL OXYGEN 

NIT = 0 
CALL PAGE 
PRINT 204, NAME 

204 FORMAT ( IOX# A8//50H NIT Z CON DELCON F(CHI) I/K 

1 FL ) / ) 

DO 12 I - 1 *NEL 

12 PRINT 206-* NIT, NZ(I), CON ( I ) 

206 FORMAT (2I5.4F10.6) 

16 TOTAL * 0.0 

DO 17 I = 1.NELP2 

17 TOTAL = TOTAL + CON(I) 

I F ( IOX ) 18,20 

18 PRINT 206* NIT* NZ(NELPl), CON(NELPl) 

I F ( IOX. GT. 0)20, 19 

19 PRINT 206, NIT* NZ(NELP2)* C0N(NELP2) 

20 PRINT 208* NIT, NTOTAL* TOTAL 
208 FORMAT ( I5*A5*F1C.6> 

LFtNIT.LT. IMIN)26*27 
27 I F ( N I T • EQ • I MAX ) 8 3 * 2 9 

2 9 DO 31 I * 1 * NVAR 

IF ( CON ( I ) .GT 0 .01, AND. DELCON! I ) .GT .DM I N ) 26 * 31 
31 CONTINUE 
GO TO 85 
26 PRINT 400 
400 FORMAT ( ) 

NIT = NIT + 1 _ - 
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DO 25 1 = 1»NVAR 
FCHI = ABSORB ( I) 

FL = FLUOR ( I ) 

CON2 = OCON( U/(FCHI*( 1*0 + FL ) ) 

DELCON ( I ) = CON2 - CON ( I ) 

IF { 0 ( I ) ) 21 *24 

21 CON(NELPl) = CON ( NELP1 ) + DELCON ( I ) *0 (I ) 

IF { SI ( I ) ) 22 ,24 - - — 

22 CON ( NELP2 ) = CON < NELP2 ) + DELCON ( I ) *S I C I ) 

2 A CON i I )- = COM2 - 

25 PRINT 206* NIT* NZ ( I ) * C0N2 * DELCON(I), FCHI, FL 
GO TO .16 

85 NS = NS + 1 

IF ( NS.LE..500 }86.*S0 
60 PRINT 260 

! 6 0 FORMAT ( / //57H.-DIMENSI0N- .EX C-EED ED - ^REMAl-NXN.q .P.O..I-NIS . .C AN NO T. BE CORRE 
1CTED ) 

62 READ 13*110) NAME 
.10 FORMAT t A8 } 

IF(E0F,3) 75*62 

86 NAME2 ( NS ) = NAME 

TOTLNS) = TOTAL . 

DO 40 J = 1*NELP2 

AO Y { NS.* J ) - CONLJ1 . _ - 

GO TO 8 

7 5 CALJ PR NT PLOT tNP* NELP2 ) . 

I F { I OX ) 50 * 1 

50 DO Z6 ) .= LUNEI ... .. . 

XX = 1.0 + 0 ( J } + SI ( J) 

. -DO - 16— L =~1*.NS 

76 Y t I * J ) = Y ( I ,J )*XX 

CALL PRNTPLOT ,1NQM*NEL) - - 

GO TO 1 

END-... - 


SUBROUTINE .CONSTANT 

TYPE REAL LAMBDA, MU 

.. .COMMON/ D AXAZ A1 32- U. LAMB DA4-9-2-)-*- EDGF 1 9 2 » 7 U Cl 91 » A) , EXPt-92) * WAVL * 

1 S INTHETA, AWK* BWK * CWK * CSIG* XSIG DATA 

DATA. (A = 1. 00S,.4. 003 *6.. 9 39, 9. 0.12 *10-.Aiaa2. Oil *14.007,15*999, A 

1 19.000,20.183*22.990,24*310,26.980,28.090*30.974,32.064,35.453, A 

2 39. 948 *39. 102 *40. 080, 44 *9.60, 47..9CKU5Q-.-94Q ,.5-2-. *54.94,55.85,58.93 A 

3 ,5 8.71 ,63.54,65.37,69.72*72.59,74.92*78 .96,79.909,83.8*8 5.47, A 

4 8.7 . 62., 8 8..9.03.^ai « 2 2.32.9-l-»-SL5-.9,4 , 98 . , 1 1 02 . 9 0 5 »AQ6 .4.107.87, A 

5 112.4,114.82,118.69, 121. 75 * 127.6, 126.9, 131 .3 , 132.905, 1 37. 34 ? A 

6 13 8. 91, ,140.1 2., 140. 91., 14.4^24 ,147. *J.5jQ...3-5-»-152.. *157.25,158.9 2, A 

7 162*5, 164*93, 167. 26*168.93,1 73 .04 ,174.9 7 *178. 49* 180.95, 183.85* A 

8 186*2* 190.23.1 92. 2, 195.09^196. 97, 2D 0^.5^, 2.0.4.37 ,207 . 19, 208. 98) 4 

DATA { A ( 84 ) = 210. ,210., 222**223., 226. ,227**232. 04, 231. *239. 03) A84 

DA T A.t LAMB DA L51 = 6 7.6 s 4 4. Z-* 31 -*6-*^.3^A2-,-ia*12-, 14. 61 * 11 .91 *9.89 *8.34 , LAMBDA 

1 7.126*6.157*5.373*4.728*4.192*3.742*3.359,3.031 *2.749,2.504, LAMBDA 

2 2.*29, 2.102, 1*937, 1.79,1.658, 1.541,1. 436 ,1.341 ,1*255,1.176*1.105 LAMBDA 

3 ,1.04,7.817) LAMBDA 

DATA ( LAMBDA ( 37 ) = 7 . 318 , 6. 863 , 6 .449 ,6 .07^5. 724 , 5 .407 , 5 . 1 1 5 ,4 . 846 , -LMDA3 7 

1 4.597*4.368 ,4. 1 54 ,3.956*3.772*3. 6 *3.^39 *3 .289 *3. 149,3 . C 17 *2. 892 LMDA37 

2 *2.776*2. 6-66 *2 ♦ 56l*-2 .443-^2-* 37 * 2-.2- 82- -*2-^2--*2 .121 *2 .047 * 1 . 977 ♦ 1*909 LMDA37 

3 *1.845*1.784, 1.727*1. 672*1. 62*1. 57, 1.52 2 *1.4 76*1. A 33, 1.391 * LMDA37 

4 1.35 1*1.313 *1.2 76* 1.313 , 1-.20-7*-!-*-! 75^-4^-909 *6 ( 0. ) *3.94] *3.827* LMDA37 

5 3.716) LMDA 37 

DATA. .( EXP ( 3 ) .. =5 -.2* 8 8 * 2.*.8.6 ,.2 •.&5j2+.8.4j-2.».83.iZ*32 *2.8 1 . 9 . 2 .80*2.79,2.79* EXP 
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1 2. 78, 2. 77, 2. 77*2. 76 *2. 76 *2. 75*2. 75 *2. 74 >2. 74 *3(2. 73), 2(2. 72), EXP 

2 3(2.71) ,3(2 .-71*3 (2.69-1 ,£>(2.-73.) *.5X2-^2-).*.4(2V714-* 342 .7 ) *3 ( 2 * 69 ) , EXP 

3 3(2.68) ,3(2.67) ,2(2.66) *3(2.65) ,2(2.64 ) ,2 (2.63) *2(2.62) *2(2.61 ) EXP 

data (EDGE = 9(15.)* 14.3*11.569,9.512*7.951*6.745,5.787,5.018* EDGE 

1 4.379*3.8 71 *3.437*3 .07 ,2 . 757 *2^497.2^.269 *2.07.1*896 ,1. 743 * EDGE 

2 1*608, 1.488 *1.38 *1.283 *1.1 96* 1.1 17, 1.045,. 98 *.92 *6.46 *5. 998* EDGE 

3 5.583,5.232 ,4. 8 67 *4. 5 8 1-*4. 2 98*4.06-^3*8-3 *-3 .626,3 . 428 ',3. 254,3. 0S5EDGE 

4 ,2.926,2.777,2.639,2.511,2.389,2.274,2. 167,2.068,1.973,1.889* EDGE 

5 1.811 * 1.735, 1-^65, 1.599-,-W-536*4*477,W4^4*l. 36-5, 3r.317*l. 268, EDGE 

6 1*222,1.182,1.14,1.1,1.061,1.025 , .99, .956 » • 923 , • 893 , • 863 , • 835 , EDGE 

7 .808 , .782 ,3.099 ,6 (U. ) ,2.401 ,2*31-3*2.235 ) EDGE 

DATA ( EDGE ( 104 ) = 16(15.), 14.24, 11.27, 10.3, EDGE104 

1 9. 54, 8.73-, .8^107,-7 .-506*6. 97, 7. 21 ,6.643, EDG5104 

1 6.172,5.755,5.378,5.026*4.71 8 *4. 436 *4 . 1 8 *3 .942 * 3 . 724 • 3 .53 4 ^ .32 EDGE 104 

2 6*3. 147*2*982 *2^33 *^687^2.5 53 *2*429*2^344-* 2. 204*2. 103*2. Oil, - FDGE104 

3 1.924,1.843,1.767,1.703,1.626,1.561,1.501,1.438,1.39,1.338, EDGE104 

4 1.288* 1.243,1. 199,1.155, 1.1 14,1.075,1.037,1.001 ,.967,. 934, .903* EDGE 104 

5 • 872, .843, .815, 3. 355, 6(0. ) ,2.577,2.479,2.394) EDGE1C4 

DATA (EDGE ( 214 ) = 11.87) 

DATA ( EDGE (215) = 10.9* 9.94 ,9. 124 , 8 .416 ,7. 8 ,7.43 ,6.89*6. 387, 5.961 , EDGE215 

1 5.538,5.223 ,4.913,4.632*4.369,4.13,3.908,3.698,3.504,3.324, EDGE215 

2 3.156, 3., 2.855, 2.719, 2.592, 2.474, 2.363, 2.258, 2.164, 2.077, 1*9 Q 5, EDGE215 

3 1.918*1.845,1.775,1.71 ,1.649, 1,579* 1.53b, 1.482,1. *33, 1.386, ECCE215 

4 1.341 *1.297,1.255,1.216,1.177,1.14,1.106 ,1.072, 1 .04 , 1 . 008 * • 979 , ED6E215 

5 .95*3.902,6(0.) ,3.08,2.98,2.882) EDGE215 

DATA ( EDGE ( 308 ) = 10 .2 * 9. 39 * 8 • 67 * 8 . *20 ( 10 • ) , 9 • 557 * 9.042 * 8 .614 * EDGE308 

1 8.188,7.831,7,513,7.178*6.856,6.566*6.292* 6.088,5.82,5.581, EDGE308 

2 5.366,5.161 ,4. 972 ,4 • 753 ♦ 4. 569 *4 • 405 ,4 . 2 31 064 ,3. 91 5 , 3. 774 , EDGE 308 

3 3 .62, 3 .4.82 *3.349 *3 .2 18,4* 808,6(0* ) ,3.745,3.614,3.4 96 ) EDGE 308 

DATA ( EDGE (401) = 3(lu.)), (EDGE(424) = 1 0 . 8 , 10 . 3 * 9 . 6^2 * 9 .2 68 , EDGE401 

1 8.773,8.376,8.023, 7.642,7.322,7.011,6.715,6.453,6.179,5.931, EDGE40 1 

2 5.686,5.475,5.211,5.015*4.823,4.625,4.443,4.273,4.114,3.939, EDGE401 

3 3.779*3.632*3. 484*10. *6(0. ). *3( 1CUJ_1_ _ „ EDGE4C3 

DATA ( EDGE ( 6 1 1 ) = 6 ( 10. ) , 9 • 957 , 9. 528 *9. 1 15 * 8 • 773 ,8 .433 ,8 • 086 * 7. 777 EDGE61 1 

1 ,7. 427 *7. 128^6. 871, 6.594, 6. 357, 6. 094,5. 8b *5. 629 *5.4 13* 5.2* EDGE611 

2 4.994), ( EDGE (518) = 1 0 .4 ,9 . 957 *97^99 *9 . 1 15 , 8 . 705 * b . 3 3 1 , 7 . 99 7, EDGE 611 

3 7. 685, 7. 383, 7. 128, 6. 834 *6. 559* 6. 331, 6.1 12 *5. 86 1,5. 637* 5.45 1 . EDGE 61 1 

4 5.253*5.06*4.873*4.71,4.522*4.349*4.201 *4.042) EDGE 511 

DATA ( C ( 3 ) = . 135-*. 35*. 74^1^5-, 2. 24- *-4^-9 +4- W 7 .9.05*11 .7-5*14.87* C3 

1 18. 5*22. 5, 27. *31. 7*36. 9*42. 5*48. 4*55. 1*62. 1*69. 8*78. ,86. 7* C3 

2 95.8,105.5, 115.9, 126.8,138., 149.8.-162.2 *175.4*189.4,20 5. ,31.25* C3 

3 33.9,36.5,39.3,42.3*45.5,48.5*52.1*55.6*59.3*63.1 *67.3*71. ,75.5 C3 

4 ,79.6, 84.2 *88. 5 *94. *98.3 *103.4 *109..* 1 14.4 *120. *125.5 *132. » C3 

5 138.5*143.5 *150. ,157. *163.5*170.5*176.5 *184. , 192 . 5 » 199 . * 206 . * C3 

6 214.,222.*231^3 239.,247.-*256. ,2 6-3^-*-2X2^-*2-81. * 289. ,298* *79.94 * C3 

7 6(0. ) *109.4*114.5*120.7 ) C2 

DATA ( C ( 1 04 ) - *89 * 1 .18 *1 .54*1 .96*2.43 *2*98*3.62 *4.31*5. 1 *5.99* C104 

1 7. *8. 02*9. 18, 10. 45 *11. 75 *13. 2b* 14. 8, 16. 45, 18. 2b* 20. 2*22. 15, Cl 04 

2 24. 25 * 26. 4* 28. 8 *26. 7 *28.8 *31. 2 *33*-5 *36. 15 *38 .9 *41. 45 *44. 5 *47. 5* C1C4 

3 50. 7,53.95, 57. 5 *60. 7 *64. 5, 68., 71.95, 75. 65 *80. 35 *8*. ,88. A ,93. 15* C3 04 

4 97*8 *102.5* 107-.3* 112 .8*118.5* 122.5* 128.. *134. * 140. * 146. » 151. * 157 Cl 04 

5 . *165. *170., 176. *183. *190 **197. *204. *211. *219. ,225**232. ,240. , 004 

6 247. ,255. ,68. 9*6(0. ) ,94*3*98. 7, 1-04. I, <C(2l4)*l5, 6) 004 

DATA t C ( 2 1 5 ) = 17.25*18.95*20.75*22.55*24.6*19.15,20.8,22.4,24. 1, C215 

1 25.95927.9*29.75,31.95,34.1*36.4,38.7*41.3*43.6*46.3,48.8*51.7* C215 

2 54.3, 57.7, 60. 3, 63. 4 *66. 9 *70. 2 *73. 6, 77. *81. *85., 88. *92. *96.3* C215 

3 100. 3, 104.5 *108.5 *113., 11 8. , 122. *126* • 1 31 • * 1 36. , 142 . * 147. * 1 52 . * C215 

4 157. *161. ,167. *172. *177. *183. ,57. 1,6(0. ) ,78.1*81.8*86.2) C215 

DATA ( C( 308 ) =13.6 *14.9* 16. 2 *17.7*5.6 *6.1 *6.62 *7. 18 *7. 78 *8. 39 *9. 06* C308 

1 9.74*10.47*11.24,12.03,12.87,13.76*14.68,15.65*16.66*17.7?* C308 

2 18.82*1 9,*.9 7 321.. 16 *.2 2 . 4 1 *23^.7 jl25.^05.32.6-»-4JE32.7^93.*.25L. 4.4 >31. 01* C308 
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3 32.63*34.31*36.06* 

4 52.44*54.79*57.21*39 

5 6(0. ) *67.5*70.7,74 
DATA ( C ( 401 ) = 4*23*4 

1 24.1,25.4*26.7*28. 

2 43.2 *45.2 *47»2 *49. 

3 6(0. ) ,31.3*32.4,33 

4 24*5,25. 8, 27 ..*28. 5 

5 42.6*44.5,46.4,48. 
DATA 2 C( 611 ) =46.3* 

1 25.7,27. ,28.3,29.6 

2 43-.6 *45 • 4 *47*3 ) * ( 

3 18. ,18.55,19.05*19 

4 24.4) 

SINTHETA = .79335 

- AWK = -.064 . - 

BWK = .034 

CWK a -1 • 03_E— £ 

WAVL = 12397.7 
CSIG = 2.39E+5 
XSIG = 1.5 

--END 


37.86*39.76*41.69*43. 
.7*62.27* 64 .-9-2 .67 .64 * 
.5) 

.69*5.43)- * (C (424j = 
1*29.6,31.1,32.6,34.3 
3 * 51 * 5 * 53. 7-+S-6.-* 58.3 , 
.6), ( C ( 5 1 8 ) = 17.9,1 
*29.8*31. 2-*22. 7 » 34 ..2 , 
3, '5 0.4* 52.6*54.7) 
17.2^-18^2- 
,30.9,32.4,33.8,35.3, 
C ( 709 ) = 1 4 .-8-, 45.3,15.- 
.6,20.1*20.45,21.35,2 


69,45.79*47.92,50.13, 

70.-55*73.6*76.64.49.3- 

19 ^ 3 , 20 . 4 , 21 . 6 * 22 . 8 , 
*35.9*37.7.39.6*41.3, 
60-^8 ,63.5,66.4*25.1 * 
8.9,20. *21. »?2.2*23.3» 
35*8,37.5,39.1*40.9, 

l+Z+2-2 ♦ 3 , 2 3 .4-* 24. 5 * 
36*9*38. 4*40.1 **1.8* 
85*-l 6.3, )6. 85, 17.45, 

1 .95,22.5,23.1 ,23.8, 


--SUBHOJJX1NE . S EXUP- ... . 

TYPE REAL LAMBDA, MU 

COMMON N PROP, _LDlL 6. ) -,—NEL - -N2_L2XLM- -SCO N (. 2 04-* 3 C F ( 2 04-* - MU BAR, 

1 0(20), IOX* NELP1* NF ( 20 ) * NZF(20,4), 2(20), R(20,4), FC(20), 

3 3-IjGMA.L 20-)-*.-WK4-20-U— MUi^_Q^2Jl)-^.-aL2il*Al^.-B.IT-2a)-*-NV.AR.* -NE-LP2 , 

3 OCON (20), CON ( 20 ) » DELC0N(20), NOM(20)* NS, NP(20), I M I N , 

4- .XMAX,.„JLGj— YKAX* UM1N- - .. 

COMMON/ DAT A/ A ( 92 } » LAMBDA (92), EDGE (92,7), C(92,8), EXP (92), WAVL 

-_1 3 INTHEXA-* AWJl,--BMiU-JCtt K. * C 3 I. G , ■ . X 3 LG . 

DIMENSION X ( 20 ) , Y ( 20 ) , CR(20), XR(20) 

-READ--10-Q-,— NP-R QBjl-IQ.,. N F I .^.-EO^.-NEJX^imN.^— 1MAX* IG, YMAX, DM IN . 

IF( EOF, 50)3*5 

. . 3 ..-S.I.QP - 

5 CALL PAGE 

R R INI— 200 

READ 102* ( NZ ( I ) , SCON(I), SCF(I), 0(1), S I ( I ) , NOM ( I) * I = 1 *NEL) 

JPJRJJXL. 202-*— QjlNZXX l * SCJQN-UJ-* SJCE.UJ.*.Q 1XX*B4 U ) -,-NCM ( I ) * 1 = 1 * NEL ) 

LF » 1 

. NYAR . = „NEk_-.NEJ.X 

DO 7 I = 1 , NEL 

LEJ 0( 1 )-lA,.7 

7 CONTINUE 

* . NELP-l- « .Jv| EL P-2- --ME L - 

GO TO 11 

, _ .. . a.;iax._= i . 

NELPl = NEL + 1 

. . NZJ.JSIELP44 =2 

DO 9 I = 1 , NEL 

. IF I.S-I-LL ) J-l Q *-9- - - -• 

9 CONTINUE 

N£LP_2 = NELR1 - 

GO TO 11 

-X-0-44EUR2-.-.--NEUP-4— +- 4 

NZ(NELP2) = 14 

— i ax. - -iox - 

11 DO 20 I = 1,NVAR 


C 303 
C308 
C 3 0 8 
^*4 0 1 
C4^1 
C 4 C 1 
C*01 
C401 
C^Ol 
C6) 1 
C611 
C61! 
CM 1 
C6i ; 


DATA 



I F ( NZ ( I ) .GE.36 J 2 0 , 1 2 

12 IF(NZU).GT. 21)15*13 

13 N1 ^ NZ ( I } + 1 
GO TO 16 

15 N1 = NZ ( I ) + 2 

16 N2 = N1 + 3 

DO 18 J r 1.NELP2 

IFtNZ ( J) .LT.N1 *OR.NZ( J) *GT*N2 ) 18*17 

17 N = N + 1 
LF = 0 

NZF ( I »N ) - J 

18 CONTINUE 

20 NF( I ) = N 

DO 40 I - 1*NELP2 

N1 - NZ ( I ) 

Z ( I ) = FLOATFtNl ) 

ENCODE ( 6 * 320 * NP ( I ) ) N1 
X = EXP { N-l ) 

IF (Nl.GT. 35)23 *21 

21 X { 2 ) = X C 3 ) - Xi4) s 2*73 
X ( 5 ) = 2*60 

GO TO 27 

23 X ( 3 ) = X’t 2 ) - X 
1F(.N1*GT*62)25*24 

24 X 1 4 ) - X ( 5 ) = 2*60 • 

GO JO 27 

25 X ( 5 ) s 2.22 
X t 4 ) = 2.60 

27 X ( 6 ) = X( 7 ) = 2*60 
X l 8 ) = 2*22 

EC ( I ) a WAVL/EDGtt N1 ) * *001 

SIGMA (I*) = CSIG/ ( EO**XSIG - EC(I)**XSIG) 

DO 35 j = 1 *NEL P2 
N2 = NZ t J ) 

Y ( J ) a LAMi3DA{N2) 

DO 30 K - 1*7 

I F ( Y ( J ) * L E • EDGE ( M 1 * K ) > 3 5 » 30 
30 CONTINUE 
K = 8 

3 5 MU ( J * I ) = C(N1»K)*YIJI*#X(K) 

IF(LF )40,38 

38 CRt I) = C(NI ) / C ( N 1 *2) 

XR ( I ) = X - X( 2 ) 

YY= (AtoK. + BWK*Z(I) + CWK*Z ( I ) #*3 ) **4 
WK( I) = YY/ ( 1 • G + YY) 

40 U(I) = EO/ECt I ) 

PRINT 500 

PRINT 208 * , NZ t I ) * I = l *NVAR ) 

PRINT 209 

DO 42 I = 1 * NELP2 

42 PRINT 210* NZtD* ( MU ( J * I ) * J=1 * NVAR ) 

CALL PAGE 
IF I LF)65*45 
45 PRINT 204 

DO 47 I = 1*NELP2 
N1 NF( I ) 

DO 46 J = 1 *N1 
N2 = NZF ( I * J ) 

46'R(I*J) -- CRt I )*Y<N2) **XR( I )' 

47 PRINT 206* I* NZ ( I ) * Nl* ( NZF ( I* J ) * R U * J ) * J= 1 » Ml ) 

PRINT 212 ' 

PRINT 214* ( NZ { I ) * Y ( I ) * EC ( I ) >SIGMA( I ) *W K( I ) >U{ I ),I = 1*NELP2) 
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GO TO 67 

-65 PRINT -24-3 — - 

PRINT 216 * lNZ(n*Y(ntEC{I)*SIGMA( I ) * I = 1 * NELP2 ) 

- 6 7.-CQNT.imiE 

PRINT 218 

.. . . .1 F_( L.F-1 J-2-*-69 . - - - 

69 DO 70 I = 1,NELP2 

. - .X.UJ =__UT-l4.*k-OGEtUiXU--«--UlI-4---^-4^Q- - - 

WK< I ) = WK ( I )**5 

-- 70 £0(44 - = P0WRF4 EC ( I ) < 0 * 3 ) - 

DO 71 I = 1 *NELP2 

-N~=- NZ(-U --- - -- -> 

M = NF ( I ) 

. . -.-Q0~71-X-^X^W - 

L = NZF(ItJ) 

K_ = .. N 2 < L ) 

71 U(I*J) = (X(L)*(R( I *J)~1«0)*WML)*ECU )*A(N)*MU< L* I) )/(X( I )*R( I * J) 

1 *EC.C L.) *AXKJ4- ... — - 

72 CONTINUE 

„„ .PR.IN T..„5'Q0_ 

100 FORMAT ( I6*6A8*I3*F3.0*4l2*2F5*0) 

J..Q.2 EO.RMAXi-5X*_L5_»AE 1 Q^6 jl54UA34 

200 F0RMAT(11X,53H STANDARD FOR ELEMENT Z OXYGEN SILIfO 

. .„lNZ£4tt M X„_. CONCENTRATION- EJJCH.1 L* 11+ FLUOR.) .FACTOR FAC 

2T0R/ ) 

-20 2_JiOJlMAJJJJ,»I^ s Fa3_%6_iF^O_^j_EX3^6^XQ^j»-5.XAA.5_)._._- . 

204 FORMAT ( / / / / 1 3X *9HNUMBER OF/HX* 11H FLUOR ESC ING/ 89H ML Z EL 

- 1EMENT. 5 EL EME NT N UMBE R A N D A BSOR PT ION J U MP-.BAI.LQ,QF-_EJU1QRESCING E 
2kEMENTS/) 

.2.0 6__EQBMAXt.l3 . juL5-*-IXQ4u41_I2_s.EB.*3-1X - 

208 FORMAT ( 10X * 37HMASS ABSORPTION COEFFICIENTS MU/RH0-//9H EmITTEI^, 

1 ±3sl&lbJL - 

209 FORMAT (9H ABSORBER) 

--■2JLQ---EQ6MA T < 1 4.»F3.*-0-» I8F6 *0 ) 

213 FORMAT {///28H Z LAMBDA EC SIGMA/) 

. .;2.12— FORMAT — (-/-/V-4-4H Z LAMBDA -EC- -S46MA- OMEGA - 11/ ) 

214 FORMAT t I 3 , F9 . 4 * F7. 3 > F9. 1 * F8.4*F8.3) 

21-6 .F.ORMA T. L13-* E9 . 4.*.EZ..a*E.9^lJ 

218 FORMAT (///37H LAMBDA = WAVELENGTH OF EMISSION LINE/ 

1 18H.. Z -.p—AJ QM I C ■■ NUMRER /'26H -R ■ ^-ABSORPTION — iU MP—iAA T 1-0/ 

2 38H EC = CRITICAL EXCITATION VOLTAGE (KV)/29H OMEGA = K FLUORES 

3jCEN C E__ YXEliDXlOH— U-J e. -EQyEC ). 

320 FORMAT ( 2 X » I 3 ) 

5.0O._FJ0RMAT.{y_/y_l — — 

END 


SUBROUTINE RAGE 

CQMM.GJNL NERQBj.. .1D-L64* _HEXj__EOa_-NiL12.Qj-t-XCON_L2lL4-5-_SCE-L2iU » MueAR , 

1 0(20), IOX. NELP1* NF ( 20 ) * NZF(20*4), Z(20)* R(20*4)* EC(20)* 

.. ? SIGMA t 70) t- WK( 20) « IH20.4). ST (20 )-* — NVAR-*. -.NELE2-*. 

3 OCON { 20 ) * CON { 20 ) * DELcON(20)* NOM(20)* NS* NP(20)* IMIN* 

..4 XMAX-*- -1jG.*-_Y-MAX-*---DMXN — - - ------ 

NPAGE = NPAGE + 1 

PRINT -2QQ-^_N2RQB+ --LD-^-EO,*. -NRAG.E 

200 FORMAT ( 12H1PR0BLEM NO* * 1 6* 3X*6A8 *5X *F2 .0 *3H- KV * IOX *4HPAGE * 1 3 / ) 
-END- 


.jSUaRQUXINE-JOX-YGEfcL. 



COMMON NPROB, I D ( 6 > » NEL * EO> NZ(20), SC0N<20), SCF(2U), MUBAR » 

1 0(20)9 IOX, NELPl, NF(20)s NZF(20»4), Z(20)» R(20,4)» EC<20)> 

2 S I GMA ( 20 ) * WK ( 2 0 ) , MU(20,20), U(20,4)» SI (20) > NVAR , NELP2 9 

3 OCON { 20 ) » CON ( 20 ) 9 DELCON<20). N0M12QJ* NS, NP<20), IMIN, 

4 IMAX 9 IG 9 YMAX, DM I N 
IF( IOX>5,2,10 

10 CON CNELPl ) = 0.0 
DO 121 = 1 9 NEL 

12 CON(NELPl) = CON(NELPl) + CON(I)*0(I) 

GO TO 2 

3 CON(NELPl) = CON ( NELP2 ) = 0.0 
DO 7 I = 1,NEL 

CON(NELPl) = CON ( NELP1 ) + C0N(I)*0(I> 

7 CON { NELP2 ) = CON(NELP2) + CON(I)*SI(I) 

2 RETURN 
END 


FUNCTION ABSORB ( N ) 

COMMON NPROB 9 I D { 6 ) * NEL 9 EO 9 NZ(20>9 SCON(20)9 SCF(20)9 MUBAR, 

1 0(20), IOX, NELP1, NF(20), NZF(20,4), Z(20), R(20,4), EC(20J, 

2 SIGMA ( 20 ) 9 WK ( 2 0 ) , MU(20,20), U(20,4), SI(20)i NVAR, NELP2 , 

3 OCON (20 ) 9 CON ( 20 ) 9 DELCON(20), NOM(20)» NS, NP(20), IMIN* 

4 IMAX, IG, YMAX 9 DMIN 

COMMON/DATA/ A ( 92) , LAMBDA (92 ) 9 EDGE (92 ,7), C(92,8), EXP (92 ) 9 WAVL* 
1 SINTHETA, AWK, BWK , CWK, CSIG, XSIG 
TYPE REAL MUBAR 9 MU 
ABAR = ZBAR = MUBAR = 0.0 
DO 10 I = 1 , NELP2 
N1 = NZ ( I ) 

MUBAR = MUBAR + MU ( N f I ) *CON ( I ) 

ABAR = ABAR + A(Nl)*CON(I) 

10 ZBAR = ZBAR + Z(I)*C0N(I) 

H = 1.2 * ABAR/ZBAR**2 

CHI = MUBAR/ ( SINTHETA*S IGMA ( N J ) -+— 1 .0 

ABSORB- ( 1 .0 + H) / ( CHI#{ } .0 + H*CHI)) 

end 


FUNCTION FLUOR ( NN ) 

TYPE REAL MU, MUBAR 

COMMON NPROB, ID(6>, NEL, EO, NZ(20), SCON(20), SCF(20), MUBAR, 

1 0(20), IOX, NELPl, NF(20), NZF(20,4), Z(20), R(20,4', EC(20), 

2 S I GMA ( 20 ) , WK ( 20 ) s MU ( ?0 ,20 ) , ±1 ( 20 .4 ) , SI(20), NVAR, NELP?, 

3 OCON (20) , C0N(20), DELcON(20), NOM(20), NS, NP(20>* I-* IN, 

4 IMAX, IG, YMAX, DMIN 

COMMON/DATA/ A ( 92 ) , LAMBDA (92), EDGE (92 ,7), C(92,8), EXP (92'. -aAVL, 

1 SINTHETA, AWK, BWK, CWK * CSIG, XSIG DATA 

N = NN 
FLUOR = 0.0 
IF ( NF (N) ) 3 ,2 

2 RETURN 

3 M = NF ( N ) 

DO 2 5 J - 1 ,M 
L = NZF(N>J) 

IF ( NZ ( L ) *GT. 35)5,7 
5 PR TNT 200, NZ(L) 

200 FORMAT (//58H FLUORESCENCE CORRECTION CANNOT BE CALCULATED FOR ELEM 
1 ENT ,12//) 

GO TO 25 



7 FMUBAR = 0.0 

DO 10 I * 1*NELP2 

10 FMUBAR = FMUBAR + CON ( I ) *MU ( L * I ) 

V “ SIGMA (L) /FMUBAR - 

V = LOGF( 1.0 + V)/V 

ULC = MUBAR/(FMUBAR*S IN-THETA ) 

ULC - LOGFtl.O + ULC ) /ULC 

FLUOR = FLUOR- +. '-CON ( L ) *U-( N* J ) * t ULC---+ VJ /-EMU BAR 
25 CONTINUE 

END - - 


SUBROUTINE PRNTPLO.T CN*.NQ.) 

COMMON NPROB* I D ( 6 ) * NEL * EO* N2{20)* SC0NC20)* SCF(20)* MUBAR* 

1 0 ( 2G.J * I OX*--. NELPLj . NFL2-XU * -NZ.F.I 2 0 *44-»-Z.l 2 0-U RL20-.4 ) » .EC ( 20 ) * 

2 SIGMAI20), WK ( 2 0 ) * MU(20»20)* U(20*4)* SI(20)» NVAR* NELP2 * 

3 OCON { 20 ) > CON (20 ) 9 DELCON(20), N0M(20-)* NS* NP(2C)* I MIN* 

4 IMAX* IG* YMAX* DMIN 

COMMON/PP/ NAME2 ( 50G ) * Y1500»20), TOTI500) 

DIMENSION I PREF ( 9 ) * X(500)* N(20) 

DATA (IPREF = A*3.»5 *0 * l-*2-*-6-*7 • 9 ) 

IF(NO*LE. 10)3*1 

1 N1 = N3 = 1 - . -- 

N2 = 10 

GO, TQ, 4 - - - 

3 N1 = 1 

N2 * NO 

4 CAL-L PAGE , 

PRINT 220 * ( N( I ) * I =N 1 *N2 ) 

220 FORMAT ( 25X * 56HCONCENTRAT I ONS CORRECTED FOR ABSORPTI-ON AND FLUORE 
1ENCE//15X*5HT0TAL » 10( 5X,A5 ) ) 

DO— 1.0 I = 1 * NS . . - - 

PRINT 256 >NAME2( I ) » TOT ( I } * ( Y ( I * J ) * J = N1 * N2 ) 

256 FORMAT ( 2X » A8 * 1 1 ( 2P F 10.4 ) ) . - 

I F ( XMOD.F ( I *5) ) 10*5 

5 IFtXMODFI I *45Ua*7 - - - - - 

7 CALL PAGE 

PR I ALT- 22-Q-* . ( N-(4C-1-*4C=N 1»N2 ) - 

GO TO 10 

8 PR-IN.T 400 ----- - - - - 

400 FORMAT ( ) 

10 CONTINUE 

IF ( N3 ) 1 1 * 14 

1.1 N1 =--11 - - 

N2 = NO 

N3 =. 0 - - 

GO TO 4 

14 IF ( I G ) 15 * 90 
90 RETURN 

15 1FI N. NE. NP ) 3Q.J.16 - - 

16 IF ( IOX ) 90 *3C 
30 YLEN . = 8.6614 

YYMAX = YMAX*1.1 
XLEN = FLOAT F t_(NS+7 ) / 8 ) 

PX1 = XLEN +1.0 

PX2 = XLEN + 1*5 — - - 

CALL CCPLOT(XLEN+3. 0*0*0) 

DY = YYMAX/YLEN 
33 DO 35 I = 1 * NS 
35 X( I ) = 1 



NX = XMI NOF ( NO * 9 ) 

P-Yl = ...8924-- 

PY2 = .7874 

DO ,40 = l.tWL 

CALL P2LINE ( X , Y ( X * J ) * NS , 1 * 1 * I PREF ( J ) * .07 * 0 .0 * 8 . 0 * 0.0 *DY ) 

CALL CCSYMBOL— (PX1 *PYl-*..2l_* I-P-RE-F-C J ) *O.Q-*.-l ) 

CALL CCSYMBOL (PX2*PY2*.2l»N(J) >0*0*5) 

PY1 = PY1 + .5 -- 

PY2 = PY2 + • 5 

40 CONTINUE- - 

CALL P2LINE ( X * TOT ,NS * 1 * 1 * 1 1 * .07 *0 • 0* 8.0 * 0 .0 *DY ) 

CALL CCSYMBOL (PXl-.PYl 

CALL CCSYMBOL (PX2*PY2* .35*5HTOTAL*0.0* 5 ) 

CALL CCAXIS ( 0-.-. 0 ..* I D.» ~48 . XLEN_*-Q._G *-0 .0 » 8 * 0- -.8.0) 

CALL CCAXIS (XLEN.0..14H CONCENTRATION , 14 * YLEN , 90 . . 0 .0 * PY, -] 2 . 
- CALL CCAXIS {0-..YLEN- » IiD»-4 8»XLEN*0.0*X).Q»A.a -.=-,8*0) 

CALL CCAXIS CO. *7. 874* I D . -48 , XLEN , 0 .0 ,0 .0 * 8 . 0 .“8.0 
CALL CCAXIS { 0 . * 0 • * 14H CONCENTRAT ION* 14 * YLEN * 90. . C. 0 * DY * 12 . 7 ) 
DO 50 I = 1 *NS 
IF ( NAME2 ( I J.EQ.lH )50,45 
45 PX = X ( I ) *• 125 

CALL CCSYMBOL {PX*YLEN*.l4* NAME 2 (I) *45.0*8) 

50 CONTINUE 

CALL CCPLOT (0*0*-3) 

END 
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APPENDIX 2: INPUT SPECIFICATIONS 


A. Data Cards for Program EMX 


Columns 


Description 


IDENTIFICATION 

CARD: 


One at the beginning of each 
problem. 


1-6 NPROB (problem number) . Up to 

six digits (all numerals) , right 
justified*, no decimal point. 


8-54 ID - alphanumeric information 

(letters, numbers, symbols) for 
problem identification. 


NOTE: NPROB and ID will appear as a heading on 

each page of output and on any graphs associated 
with this problem. If NPROB — 0, the program 
will terminate. 


NEL (total number of elements 
analyzed in this problem). 

NEL _< 9, 

NDT (number of readings in each 
deadtime data point). NDT 99. 
Right justified, no decimal point. 
If no deadtime data, NDT = 0. 

NS (number of readings m each 
sample peak, sample background 
or standard background data 
point). NS <_ 99. Right justi- 
fied, no decimal point. If the 
problem involves only the cal- 
culation of dead time, NS = 0. 

NOTE: If readings are not grouped according to 

NS and NDT specifications, program will print an 
error message and move on the next problem, if 
any, without making any calculations; or, if the 
total number of consecutive equivalent readings 
happens to be evenly divisible by NS (or NDT for 
deadtime readings) , it will group data incorrect- 
ly. Standard peak data points may be made up of 
a variable, unspecified number of readings 100. 

If one data point includes both standard peak and 
standard background data, it will be handled as 
"standard peak". 

62 MBKGD (method of taking back- 

ground readings). If MBKGD - 0, 
the program expects to find a 
pair of background data points 
for each sample or standard peak 
data point. If MBKGD = 1, one 
pair of background data points is 
averaged and the result subtracted 
from all of the appropriate peaks 
until the time when another pair 
appears or until the end of the 
data set. If no background points 
are given, no background will be 
subtracted from the points in that 
data set. There will be no indica- 
tion of this in the final output. 


*"Right- justified" means that the rightmost digit 
of a number is in the farthest right column allow- 
ed for the number; any blank columns are on the 
left of the number. 


56 


57-58 


59-60 


Columns Descrip tion 

CAUTION: Background data points are assumed to 

occur in pairs (not necessarily adjacent), which 
are averaged to find true background. Reading 
down each spectrometer column, the first two back- 
ground points (wherever they may be) are assumed 
to make up the first pair, and so on. If there is 
an odd number of background points, the last one 
will be used alone as a true background. This 
means that if one member of the first pair is 
omitted, all further points will be grouped in- 
correctly. 


In BKGD Method 0, the omission of a background or 
peak data point will cause all succeeding background 
points to be subtracted from the wrong peaks. 

62-65 SLIMIT (limit of scatter allowed 

within a standard peak data point). 
Either right-justified integer or 
the decimal point must be included. 
If in any standard peak data point 

a > SLIMIT* 

(where o = root mean square, N = 
number of counts) 


the computer will print an error 
message and move on to data for the 
next problem, if any. It is sug- 
gested that SLIMIT = 2 to 3. 

67 IG. If IG = 1, computer will make 

a plot of the concentration of 
each element or oxide and total 
concentration of all elements 
analyzed at each point on the 
sample. If IG = 0 or blank, no 
plot will be made. 

68-70 XMAX — (the value of the current 

at the right-hand edge of the 
dead time curve plot) . Whenever 
deadtime is calculated, a plot 
is made of counts /current as a 
function of current, and a value 
must be given for XMAX. For ease 
in reading values from the plot, 
XMAX should be a number with only 
one significant figure. Decimal 
point required. If no dead time 
data, may be left blank. 

71-73 YMAX (maximum y- value on concen- 

tration plot, which will be in- 
creased by 10 percent by the pro- 
gram) . Decimal point required. 

If IG = 0, may be left blank. 
Example: if elements analyzed are 
expected to total 100%, set YMAX 
= 1.0. Plot will show concentra- 
tions up to 110%. 


ELEMENT CARD: One "element card" is needed for 

each element analyzed in the problem. Cards must 
be arranged so that NL's (see below) are consecu- 
tive integers from 1 to NEL. If elements analyzed 
are not numbered this way in probe output, NEL must 
be changed and dummy cards must be made for the 
missing elements. For example, if only three ele- 
ments are being analyzed but these are numbered 2, 
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Columns Description 

3, and 6 in the probe data, NEL must be set equal 
to 6 and there must be an element card (with NL 
in column 1) for Elements 1, 4 and 5 as well as for 
the elements analyzed. (This will take more computer 
time than if the elements analyzed had been numbered 
1, 2 and 3 in the first place.) 

1 NL (the element code number 

assigned to this element at the 
time probe data was taken) . Ele- 
ments may be numbered in any 
order, without regard to atomic 
number, order in which analyzed, 
etc. 

2-4 NO(NL) (atomic number of the ele- 

ment) . Right-justified, no decimal 
point. 

5-10 WTPC(NL) (concentration of ele- 

ment in its standard). NOTE: If, 
for example, standard contains 25 
percent Fe, WTPC = .25. Decimal 
point required; may appear any- 
where within the field. 

11-25 TAU(NL) (deadtime in seconds) . E- 

format; exponent of 10 must be in 
Col. 25; decimal point required. 

If TAU is to be computed by the 
program, put a zero somewhere 
within Cols. 11-22. If deadtime 
correction is to be omitted, put 
any negative number (with decimal 
point) somewhere within Cols. 11- 
22. If the TAU from some other 
element is to be used for this 
element, set TAU = NL cf the ele- 
ment from which the TAU is to be 
used. 

26-35 0(NL) (Oxygen factor). The oxide 

concentrations will be calculated 
as C0N(NL)*(1 + 0(NL)). Decimal 
point required; may be left blank 
if no oxygen. Example: to convert 
silicon concentration into concen- 
tration of SI02» 0 = 1.139. 

36-45 SI(NL) (silicon factor). If sili- 

con is not analyzed for, silicon 
concentrations can be calculated 
in same manner as oxygen concen- 
trations. Decimal point required, 
or may be left blank. 

51-55 NOM(NL) (name of oxide or other 

element combination to be used in 
headings on printout and plots. 
Example: A^O^*) 

PROBE DATA RECORDS: One card or magnetic tape re- 

cord for each reading from probe. If records are 
on magnetic tape, they must be written in BCD. 


Columns Description 

PROBE DATA RECORDS (Continued) 

No decimal 
points. If 
a column is 
left out in 
the output 
from probe, 
it must con- 
tain zeros 
or blank 
spaces. 

31-36 KODE (a set of six code digits to 

identify the element and type of 
reading. NL = identifying number 
of the element being analyzed for, 
and may be any arbitrary number 
from 1 to 9. If NL = 0, that read- 
ing for that scanner will be dis- 
regarded. KD = type of reading. 

If KD = 0, sample peak 

1, sample background 

2, standard peak 

3, standard background 

4, deadtime reading 

Col. 31 KD) 

32 mi J for Scanner 1 

33 KD) _ 0 

34 NLJ f ° r Scanner 2 

35 KD3 

36 ml I for Scanner 3 

37-44 NAME (alphanumeric information 

for identification of a data point) 
May be left blank. 

BREAK CARD: Data sets must be separated by a "break" 

card, which contains zeros in Columns 31-36. Any 
data in other columns will be ignored. 

After the last data card for each problem, there 
must be a card with 999999 in columns 31-36. 

After the 999999 card of the last problem for the 
computer run, there must be one blank card. (Or if 
data records are on tape, blank card must be after 
last ELEMENT CARD of last problem.) 


B. Data Cards for Program EMX2 

Columns Description 

IDENTIFICATION CARD: One at the beginning of each 
problem. 

1-6 NPROB (problem number). Up to six 

digits (all numerals) , right justi- 
fied, no decimal point. 


1-6 Time (in hundredths 

of sec.) 

7-12 Current 

13-18 Scanner 1 counts 

19-24 Scanner 2 counts 

25-30 Scanner 3 counts 


7-54 ID - alphanumeric information 

(letters, numbers, symbols) for 
problem identification. 
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_Columns Description 


Columns 


Description 


NOTE: NPROB and ID will appear as a heading on 

each page of output and on any graphs associated 
with this problem. If DATA RECORDS (see below) 
are on punched cards, only one problem may be run 
at one time. If DATA RECORDS are on magnetic tape, 
any number of problems may be run at once, as long 
as DATA RECORDS are all on one tape. In this case, 
data cards should be followed by one blank card. 

55-57 NEL - number of elements for each 

sample to be read from DATA RECORDS. 
NEL _< 20. 

58-60 E0 - operating voltage in KV. Right 

justified, no decimal point. 

61-62 NFIX - number of fixed elements. 

Concentrations for these elements 
will be read from data records 
and used in corrections but will 
not themselves be corrected. All 
"fixed elements" must appear after 
the other elements on ELEMENT CARDS 
and DATA RECORDS. Right justified, 
no decimal. 

63-64 IMIN - minimum number of itera- 

tions. Right justified, no decimal. 

65-66 IMAX - maximum number of itera- 

tions. Right justified, no decimal. 

68 IG. IG = 1 if a graph is to be 

plotted, IG = 0 if no graph de- 
sired. 

69-73 Maximum Y- value on graph (which 

will be increased by 10 percent by 
the program). If IG = 0, may be 
left blank. Decimal point requir- 
ed. 

74-78 DMIN. The calculation will be 

iterated at least IMIN and no more 
than IMAX times. Between these 
limits, iteration will stop when 
concentrations of no elements pres- 
ent in amounts greater than one 
percent change more than DMIN from 
one iteration to the next. 

ELEMENT CARD: One needed for each element in the 
problem. Cards must be arranged in same order as 
element concentrations appear on DATA RECORDS. No 
element card is needed for oxygen, the concentra- 
tion of which will be calculated on the basis of 
the other element concentrations. (Silicon can be 
calculated in the same manner; then no element 
card or probe data would appear for silicon^ 

6-10 NZ(I) - atomic number of element 

i. Right justified, no decimal. 


21-30 SCF(I) - correction factor for the 

standard of the ith element: f(chi) 
. (1 + Ifi/l-j). If SCF = 0, all SCF 
for this problem must = 0. In this 
case correction factors are calcu- 
lated for the standard, and new 
ELEMENT CARDS are punched for later 
use (see Sec. IIIC) . Decimal point 
required. 

31-40 0(1) - oxygen factor: the oxygen 

concentration in the sample will 
be calculated as ECON-jOi;!* 1,NEL. 
Oxide concentrations will be cal- 
culated as C0N^(1 + 0^) . Decimal 
point required, or may be left 
blank if no oxygen. 

41-50 SI(I) - silicon factor. If silicon 

is not analyzed for, silicon con- 
centrations can be calculated in 
same manner as oxygen concentra- 
tions. Decimal point required, or 
may be left blank. 


56-60 NOM(I) - name of oxide or other 

element combination to be used m 
headings on printout and plots. 
Example: A^O^ 

PROBE DATA RECORDS: any number up to 500. May be 

on punched cards or on magnetic tape on Logical 
Unit 3. 


9-16 

17-24 

25-32 

33-40 

41-48 

49-56 

57-64 

65-72 

73-80 


NAME - alphanumeric information 
to identify the sample or sample 
point. May be left blank. 

Observed concentrations (I^/I(A)) 
of elements in sample. Must appear 
in same order as element cards, 
with fixed elements last. E8.0 
format - for example, 33.3 percent 
would appear as 3330E-04. 

5 ppm would appear as 0005E-06. 


If the problem includes more than 9 elements, two 
records (cards) may be included for each sample or 
sample point. The concentration of the 10th ele- 
ment would appear in columns 9-16 of the second 
card; the 11th, in columns 17-24 etc. 


C. Special Tapes 

To make EMX and EMX2 as flexible as possible 
in input and output operations, the use of several 
special tapes is built into the programs. These 
tapes must be provided, bypassed, or equivalenced 
to some standard 1/0 unit; they cannot be ignored. 


11-20 SC0N(I) - concentration of element 

I in its standard. Decimal point 
required. 


Data from the microprobe is read into Program 
EMX from Tape 2. If probe data is actually to be 
on punched cards. Tape 2 must be equivalenced on 
the appropriate control card to the standard input 
unit. 
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Results from Program EMX that are to be used 
as input for Program EMX2 will be written on Tape 
3. If the user prefers to receive this output in 
the form of punched cards, Tape 3 may be equival 
enced to the standard punch unit or cards may be 
punched later from the tape. If both programs 
are being run together as a multiple-execution 
job, this intermediate output could most economi- 
cally be on a standard scratch tape; then Tape 3 
would be equivalenced to the standard scratch tape. 
This tape must be rewound between the execution of 
EMX and EMX2. If EMX is being run alone and ho 
such output is desired because no further correc- 
tions are to be made to the data, the program user 
must specify that Tape 3 be bypassed. 

If input to Program EMX2 is to be entirely from 
punched cards, Tape 3 must be equivalenced to the 
standard input unit. 

A special' tape must be provided for plotting 
information when plots are desired. 



EMX PROGRAM 
( Finds first approximation 
using data directly from 
probe) 



IG is a flag to bell 
if there are any plots 
to be made. 

Tape on Unit 3 (if not 
bypassed) will contain 
first approximations 
to be punched on cards 
for future use or pro- 
cessed immediately. 

©For explanation of 
indices of the do- loops 
involved, see charts 
for subroutines EMXINPUT 
and BKGD. 



Assign relative times . Each data set 
begins at T=0. For each data point, 
TIME = seconds elapsed from start of 
data set to midpoint of readings for 
that data point. © 
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LAST = index of last line 
of data for this problem. 

NN = element number 


Note: if TAU<0, it 

means don't correct for 
deadtime. I f TAU = 0, 
find correct TAU using 
data. If TAU>1 , use 
TAU from element tau. 


©Once a 
is found. 


sample no. NN 
it is assumed 


that all of that element 


will be on the same 


scanner . 




Ln 

U) 







Oi 

-p~ 


For oil data sets on all scanners: 
Find first standard peak in set for 
each element. Concentration ~ 

c °N = ^ IS - WTpc (NL ) 


(Dl_ook for first standard 
peak (II). If KO = 2, 
data point is a standard 
peak . 

Subroutine ORIFT interpo- 
lates between II and 12 
to correct for drift. 

If IFLAG = 1 after return 
from ORIFT subroutine, 
no drift correction was 
made between II and 12. 

Now we must check for drift 
between the original II and 
the next standard. 

UTPC(NL) = concentration of 
element NL in its standard. 


For all lines 3 all 
scanners, change ERR 
to relative error; ie., * 
ERR= ERR/COUNTS 


X Call \ 
/ BKGD(l) \ 
/ BKGD(3) \ 

f to subtract bach- j 
1 ground from sam-J 
\ pies (KD = I) and J 
\ standard / 
"^\{KD=3K/ 

















Print I record for each data point \ 
for each scanner: line index of \ 
data point, element no., sample J 
name, xurrent, counts/ sec, CON, E J 


/AreV 
^xThere the^s. 
/some no of samples''-' 
.. (or points on sam- 
\ ple)foreocti / 
^v^lement^^ 


Print error 
message 



Ln 

Ln 






















Cn 

"4 










Ln 

00 









SUBROUTINE DEADTIME (N) 



8(1) and 0(2) are coeffici- 
ents of x where f(x) = B(l) 
+ B(2)X. Calculated by LSQ 
subrouti ne. 






Subroutine DRIFT (I IFLAG) 


ENTRY 

DRIFT 




KI-K2 ore limits within 
which to correct 
for drift 

AT = TIME(I) -T 


O' 

o 








SUBROUTINE BKGD(KODE) 


J = no. of scanner 

Subroutine is done twice, 
once for K0DE=3 (standard 
bkgd) and once for K0DE=1 
(sample bkgd) . 

Ml and M2 are indices of 
first and last lines of a 
data set 


(NO BACKGROUND FOR THIS OATA SET ON THIS SCANNER ) 


Jt is first line to look 
on for a sample from which 
to subtract 8. 

O^Look for an appropriate 
background reading. 

11 = line of first of a 
background pair. 

NL 1 = element to which II 
appl i es. 

3 - true background 

2 

EBSQ. ~ (error in B) 

12 = line of second of a 
oackground pair. 

j2 is last line on which 
to look for a sample from 
which to subtract B 

Samples to which B ap- 
plies must be before any 
more background readings, 
so find 3rd background, 
i f any. 



cK 


H* 





HBKGO = method by which 
bkgd readings were ta|ten. 

®w« have already found 
the first of the next back- 
ground pair, so go back and 
find the second, if any. 


No more bkgd 
YES readings in this 
data set, so: 

J 2= M2 


KD{I,J)=KODE 

^AND-NL(I,J)=NU 

I J 2 = 1-1 | 

i - 

COUNTS{I 1,J)+C0UNTS{I2,J) 

- - 2 

ERR(II,J) 2 +ERn(I2 t J)* 


Print error messoge: 

NO SAMPLE DATA FOR 
BKGO, LINES ~~ * 


' a iwitV 

>|Q ;i t< ))to corrtct with-^ 


JVtBKGD = 0? 


Subtract 8 and compute ERR 
for sample (I, J) only 
Reset J! =1 + 1 


Subtract B and compute ERR 
for aM samples of elementNLI 
from Ji = J2 + 2 
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CONTINUE 












ON 

u> 







PROGRAM EMX2 


Store SCON in CON 



Oo for all elements, I to 

nel: 

FCHI= ABSORB (I) 

FL = FLUOR (I) , 
SCF(I)s FCHI (l.O + FL) 


f Print results of SCF \ 
calculation and punch I 
\ new "ELEMENT CARDS" / 
\ for future use / 


START 



have 

un SCF for this 

-''problem been calculated ’ 
previously 


Set NS = 0 


( Read from TAPE 3, NAME, 
OCON (I), I « 1, NEL 


OCON(I) = original concentration 
(from EMX Program) of element I 
in sample 

C ON ( I ) = concentration of element 
I in sample 

SCON(l) = concen tration of ele- 
ment I in its standard 
SCF(I) is correction factor 
[ F ‘*> ( i+r f ,/r)j 
for element I standard 

NP is BCD array of atomic nunt- 
oers to be printed as column 
headings. NOM is BCD array of 
compound names to be printed as 
column headings. 

ABSORB and FLUOR are function 
subroutines. 

NAME is code name of sample or 
point on sample; may be blank. 


NS - number of samples or points 







Set NIT = 0 
Calculate TOTALof 
all concentrations, 

I to NELP2 


Change all element concen- 
trations: 

Y{I,J)=Y(I,J)*aO+0(J) 

+ SI(J) 


Print concentrations and 
TOTAL 


XCaliX 

/prntploA 

V(NOM,NEL)/ 


NELP2 = number of elements given 
plus oxygen and silicon, if pres- 


^-^DELCON > 

OMIN for any elements 
present in concentrate* 
Ions > |% 


Store NAME in 
NAME 2 (NS), yES 
TOTAL in TOT(NS), * — < 
and CON( J) in Y 
(NS,J),J=1,NELP2 


I 


Print error 
message: 

"DIMENSION EX- 
CEEDED" 


Skip over rest of 
data for this [ 

problem ! 


NIT = NIT + I 


Do for all variable concentrations, 

[ = 1, nvar: 

FCHI = ABSORB (I) 

FL= FLUOR (I) r , 

C 0 N 2 = 0C0N(X)/JFCHr * (I.O + FLlj 
DELC0N(I)=C0N2 - CONI I) 
CON(NELP0=CON(NELPI)+DELCON( I) 
* 0 ( 1 )® 

C0N(NELP2)= C0NINELP2 + DELCON 

{I)*SI(D ® 


Print NIT, NZ, C0N2, DELCON, 1 
FCHI, FL for all elements and 
TOTAL J 


ININ and I MAX = minimum and maxi- 
mum number of iterations to be 
performed. OMIN = largest ac- 
ceptable change in concentration 
of major elements from one itera- 
tion to the next. 

® These steps done only if 
oxygen and/or silicon factors 
present. 


as 

m 








ON 

ON 










Lambda array contains the one 
wavelength (either k, L, or H) 
used for each element. 

HU = matt absorption coefficient, 
£R * ratio of Heinrich eoefffci- 
R, the abtorption junto ratio, 

(T)N0T£: after exit from this tub- 

routine, U It an expression 
involving many variables which 
Is to be used in fluorescence 
corrections. IT IS NO LONGER 

= eo/ec. 


Colculote array of obsorpfion 
jump ratios R with indices I 
■ and Jcorr«spondmg to those 
of NZF 

R(l.J)«CRU)XY(N2)XXXR(I) 
where N2 » NZF(I.J) 


f Print all Y.EC, SIGMA, and 

f Print NFUUNZFU.J), \ 

l (if LF*0)WK,U y 

l R(I t J),J=l,NF(I)),l*|,NELP2 ) 


I Calculate array U with indices I and J 
corresponding to those of NZF 


^ UiiogUt-Ut+l r I(L 2 At ^ 
where L*NZFU,J), N*NZ(l),ond K = NZ(L) © 


(RETURN) 













